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SECTION  1 


INTRODUCTION  AND  SUMMARY 


This  15-month  applied  research  program  was  aimed  at  determining  the  feasibility 
of  making  subarc-second  quality  measurements  of  earth  crustal  deviation  and 
earth  polar  axis  wander  using  ring  laser  gyro  (RLG)  based  systems. 

The  program  comprised  analyses  of  the  performance  of  two  selected  system 
configurations  and  rate  table  laser  gyro  testing.  This  final  report  describes 
the  work  done  and  the  results  obtained,  presents  conclusions  on  the  feasibility 
of  making  subarc-second  measurements  of  earth  polar  axis  (EPA) /crustal 
deviations,  and  frames  recommendations  for  further  work. 

The  two  configurations  selected  consisted  of  a  three-axis  navigator  system, 
using  three  laser  gyros,  three  accelerometers,  and  compound  gimballing;  and  a 
one-axis  system,  comprising  one  laser  gyro,  two  tiltmeters,  and  indexation  to 
point  the  gyro  alternately  east  and  west. 

The  laser  gyro  selected  for  both  the  analyses  and  the  testing  phases  of  this 
program  was  the  Honeywell  GG1389  gyro.  This  state-of-the-art  gyro  is  probably 
the  finest  in  the  world. 

Preliminary  analysis  showed  that  it  is  impossible  to  distinguish  earth  crustal 
deviation  around  the  vertical  from  EPA  wander  along  the  east-west  axis,  using 
self-contained  inertial  measurements. 

Preliminary  calculations  based  on  GG1389  gyros  with  computer-driven  random 
drift  reduction  showed  that,  in  a  kinematically  benign  environment,  combined 
crustal  deviation  and  polar  axis  east-west  wander  could  be  determined  with 
subarc-second  accuracy  over  less  than  24  hr. 


Based  on  these  preliainary  resalts,  we  defined  a  reasonable  performance  goal: 
determine  the  16-hr  average  valne  of  combined  earth  crnstal  deviation  and  EPA 
deviation,  within  0.2  arc-sec  rms,  over  a  16-hr  period,  in  a  benign 
environment . 

Using  covariance  analysis  as  needed,  a  linear  error  propagation  model,  and 
known  and  extrapolated  component  characteristics,  oar  analyses  showed  that  the 
simple  one-axis  system  should  perform  as  well  as  the  three-axis  navigator 
system,  and  that  both  coaid  come  close  to  reaching  the  performance  goal.  Rate 
table  tests  on  Honeywell  GG1389  RLGs  showed  that  this  gyro's  performance  is  now 
almost  compatible  with  the  performance  goal.  The  crucial  characteristics  are 
random  drift  and  input  axis  stability. 

The  random  drift' was  about  0 .000045°/ rt-hr,  when  operated  in  the  scale-factor 
mode  (i.e.,  continuously  rotated  without  dithering).  The  input  axis  stability 
ranged  from  0.004  to  0.04  arc-sec,  relative  to  the  gyro  body  structure,  and 
from  0.25  to  0.90  arc-sec  across  removals  and  remountings  of  the  gyros. 

The  low  random  drift  performance  could  be  retained  when  operating  in  the 
dithered  mode  by  using  a  real-time  computer-controlled  laser  gyro  random  drift 
improvement  technique,  called  CLIC,  invented  by  Honeywell.  Although  still  not 
operative  under  general  input  rate  environments,  (XlCwill  reliably  improve 
dithered  laser  gyro  random  drift  under  benign  environments,  down  to  near 
quantum  limited  performance. 

It  is  reasonable  to  suppose  that  careful  design  of  an  indexed  system  should 
improve  the  input  axis  stability,  from  the  above  0.25  to  0.90  arc-sec 
mount/ remount  value,  down  towards  the  0.04  arc-sec  intrinsic  value.  This  could 
be  achieved  by  improving  the  quality  of  the  mounting  geometry,  and/or  by 
directly  measuring  the  mounting  misalignment  and  compensating  the  system 
output . 


The  tiltmeters  used  in  the  one-axis  system  would  contribute  negligible  added 
error  to  that  from  gyro  random  drift.  Tiltmeters  are  readily  manufactured  with 
subarc-second  resolution  and  multi-hour  stability.  In  the  three-axis 
navigator,  the  accelerometers  would  make  a  larger  error  contribution  (it  is 
more  difficult  to  make  ultrasensitive  and  ultrastable  accels  than  tiltmeters), 
but  still  minor  compared  to  that  from  the  gyros. 

Considering  the  foregoing  results  and  arguments,  we  conclude  that  it  is 
feasible  to  make  subarc-second  measurements  of  combined  EPA  deviation  and  earth 
crustal  rotational  displacements,  over  less  than  one  day,  using  a  RLG-based 
system,  in  a  configuration  as  simple  as  the  one-gyro-plus- two— til tme ter  device. 

We  recommend  that  the  Air  Force  Geophysics  Laboratory  (AFGL)  design  and  build  a 
one-axis  demonstrator  unit,  as  described  herein,  to  be  tested  by  AFGL  under 
field  conditions,  as  a  prelude  to  designing  and  building  a  more  elaborate  unit 
with  all  features  needed  to  function  as  an  autonomous  transportable  field  model 


sensor. 


SECTION  2 


PRELIMINARY  STUDIES 


ASSUMPTIONS 

We  defined  the  set  of  assumptions  under  which  the  study  would  proceed,  as 
listed  below. 

o  Rigid-earth  geography:  invariant  distances  between  all  surface 

features  recognized  as  fixed  in  the  everyday  sense.  Terrains  such  as 
glaciers,  icefloes,  etc.,  excluded.  Intercontinental  drifts  considered 
negligible. 

o  Despite  the  rigidity  as  defined  above,  local  short-range  twist 

deformations  occur.  These  are  to  be  determined  to  within  subarc-second 
accuracy,  over  less  than  one  day. 

o  Gravity  vector  at  any  location  is  constant,  expressed  in  a  geographic 
reference  frame  affixed  to  the  rigid  earth. 

o  Earth  polar  spin  axis  will  deviate  from  nominal  orientation,  as  seen  in 
rigid-earth  reference  frame.  This  deviation  is  to  be  determined  with 
subarc-second  accuracy  over  less  than  a  day.  Polar  axis  wander  during 
the  measurement  period  is  negligible. 

o  Crustal  twist  and  polar  axis  deviation  to  be  determined  using  a 

substantially  self-contained  device  (e.g.,  without  need  for  star  shots, 
triangulation,  etc.). 

o  Benign  environment:  fixed  position,  zero  velocity,  and  zero  kinematic 
acceleration,  all  specific  force  is  gravitational. 


CONFIGURATIONS  AND  PRINCIPLES  OF  OPERATION 


We  selected  two  system  configurations  for  study.  One  of  these  was  a 
high-precision  three-axis  inertial  navigator,  with  three  RLGs  and  three 
accelerometers.  The  other  was  a  stripped  system,  designed  to  exploit  the 
advantages  of  operation  under  ultra-benign  kinematic  environment.  It  consists 
of  a  single  gyro,  whose  input  axis  lies  along  nominal  east-west,  and  two 
tiltmeters,  which  measure  tilt  around  the  north-south  and  the  east-west  axes, 
with  180°  indexing  around  nominal  vertical,  to  enable  bias  cancellation.  The 
principle  behind  both  systems  is  to  define  earth  crustal  twist  and  EPA 
deviation  in  terms  of  a  reference  frame  defined  by  the  local  vertical  and  the 
nominal  EPA. 

The  accels  (or  tiltmeters)  define  the  local  vertical  (unit  g-vector)  as  seen  in 
the  body  frame  of  the  device.  The  angle  sensor(s)  define  the  total  rotation 
vector  (earth  spin  +  body  rate),  also  as  seen  in  the  body  frame.  North  is 
defined  as  lying  along  the  horizontal  component  of  earth  spin,  which  is 
nominally  of  value  E+cos(LT),  where  E  is  full  earth  rate  (about  15°/hr),  and  LT 
is  latitude. 

Body  tilt  around  the  north-south  axis  through  an  angle  A(NS)  will  add  a  spin 
component  along  east-west  of  value  E*sin(LT) *A(NS) ,  comprising  the  rotated 
vertical  component  of  earth  spin.  This  will  deflect  the  total  spin  vector 
east-westwards  through  an  angle  of  (E*sin(LT) *A(NS) ) / (E*cos(LT) ) ,  * 

A(NS) *tan(LT) .  This  is,  of  course,  easily  corrected  for  by  using  the 
north-south  tilt  angle,  A(NS),  which  is  arrived  at  by  interpreting  the 
accel/ til tmeter  readings. 

Body  rate  around  the  east-west  axis  of  value  R(EW)  will  sum  vectorially  with 
the  earth  spin  vector  to  produce  a  total  spin  vector  whose  horizontal  component 
is  deviated  east/westwards  from  north,  through  an  angle  of  R(EW) /(E*cos(LT) ) . 


It  will  also  rotate  the  sensed  local  vertical  at  rate  R(EW)  around  east-west. 
This  can  be  seen  in  the  rate  of  change  of  the  accel  signals,  and  thus  east-west 
tilt  rate  can  also  be  corrected  for. 

However,  deviation  in  the  horixontal  component  of  earth  rate  cannot  be  detected 
when  using  only  the  g-vector  in  addition  to  earth  rate.  The  effect  of,  say,  a 
1  arc-sec  rotation  of  the  horizontal  part  of  earth  rate,  from  true  north  to 
north  +  1  arc-sec  east,  cannot  be  distinguished  from  the  effect  of  a  rotation 
of  the  body  on  which  the  gyro  is  mounted,  through  an  angle  of  1  arc-sec 
counterclockwise  around  the  vertical. 

This  is  the  distinguishabil ity  problem:  how  does  one  determine  EPA  deviation, 
while  using  the  EPA  as  one  ''leg"  of  the  reference  frame  in  which  the 
deviation  is  to  be  expressed?  (The  other  ''leg"  is  the  local  vertical.) 
Obviously,  the  only  detectable  kind  of  EPA  deviation  would  be  that  causing  a 
change  in  the  angle  between  the  EPA  and  vertical,  i.e.,  meridional  deviation. 
East-west  deviation  would  be  undetectable. 

BASIC  ANALYSIS  OF  SYSTEM  OPERATION 

Let  the  earth  spin  rate  vector,  including  deviation,  be  denoted: 

e  =  E  +  d, 

where : 

e  =  deviated  earth  spin  rate  vector 
E  =  nominal  earth  spin  rate  vector 

d  *  deviation  vector,  whose  direction  denotes  the  axis,  and  whose  magnitude 
denotes  the  amount  of  the  spin  rate  deviation 

Let  the  crustal  warp  be  denoted  by  a  small  rotation  vector,  a,  whose  direction 
represents  the  axis,  and  whose  magnitude  represents  the  angle  of  warp.  Let  the 


crustal  warp  rate  be  denoted  as  the  rotation  rate  vector,  p,  whose  direction 
represents  the  axis,  and  whose  Magnitude  denotes  the  aaiount  of  the  warp  rate. 

Rigid-earth  coordinates  are  defined  in  tens  of  unit  vectors  i  along  east,  j 
along  north,  and  k  along  the  upward-vertical.  The  nominal  earth-spin,  E,  is 
expressed  as 

E  =  R*(j*cos(L)  +  k*sin(L)) 

where : 


R  =  nominal  spin  rate  (about  15 .04107°/hr) 

L  =  nominal  latitude 

Similarly,  we  denote  the  warp  angle  vector,  a,  as 

a  *  i*al  +  j*a2  +  k*a3 

the  warp  rate  vector,  p,  as 

p  =  i*pl  +  j*p2  ♦  k*p3 
and  the  earth  spin  rate  deviation,  d,  as 
d  =  i»dl  +  j*d2  +  k*d3 

Further  assume  that  the  environment  is  so  benign  that,  averaged  over  the  period 
of  observation,  all  that  the  accels  sense  are  components  of  the  gravity  vector, 
seen  along  the  sensing  axes  of  the  accels.  Then,  we  can  substitute  tiltmeters 
for  the  accels.  The  tiltmeters  will  simply  (and  sufficiently)  indicate  the 
direction  only  (no  magnitude)  of  the  gravity  vector,  and  they  will  do  so  more 
accurately  than  accels  would.  Now  suppose  that  the  following  quantities  are 


determined  over  the  period  of  measurement,  as  perceived  in  the  crnstally  warped 
body  axes: 

Average  net  sensed  rotation  rate,  e 

Average  sensed  unit  gravity  vector,  g  « 

Average  rate  of  change  of  sensed  gravity  vector,  g' 

By  simple  kinematics,  potting  the  true  value  of  the  gravity  vector  equal  to  k: 

e  =  E+  d  +  p+  axE, 
g  ■  k  +  a  x  k 
g'“  P  x  k 

where  x  denotes  the  vector  cross-product  operator. 

Substituting  the  vectorial  component  representations  for  E,  d,  p,  and  a,  we 
get : 

e  =  i*(R*(a2*sin(L)  -  a3*cos(L))  +  dl  +  pi) 

+  j * (d2  +  p2  -  R*al*sin(L)) 

+  k*(d3  +  p3  +  R*al*cos(L) ) 

g  =  i*a2  -  j*al  +  k 

g'  *  i*p2  -  j*pl 

The  components  seen  along  the  three  axes  i,  j,  and  k  are,  for  the  rate  sensing 
observation: 

el  *  R*(a2*sin(L)  -  a3*cos(L))  +  dl  +  pi 
e2  -  d2  +  p2  -  R*al*sin( L) 
e3  *  d3  ♦  p3  +  R*al*cos(L) 
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for  the  tilt  observations : 


gl  -  a2 

g2  -  -si 

g3  -  1 

and  for  the  tilt-rate  observations: 

t'l  -  P2 
g'2  -  -pi 
g'3  *  zero 

Fro*  these,  we  get  the  partial  set  of  solntions: 

al  -  -g2 
*2  «  gl 
Pi  «  ~g '2 
p2  ■  g'l 

Substituting  in  the  rate-sensing  observations,  we  get: 

dl  -  R*a3*cos(L)  -  el  -  R*gl*sin(L)  +  g*2 
d2  •  e2  -  R*g2*sin( L)  -  g'l 
d3  +  p3  *  e3  +  R*g2*cos(L) 

This  shows  that  we  can  determine  al ,  a2 ,  pi,  p2,  and  d2,  but  that  we  can  only 
determine  d3  and  p3  in  the  combination  d3+p3,  and  dl  and  a3  in  the  combination 
dl-R*a3*cos(L) .  In  some  situations,  we  may  know  a  priori  that  the  vertical 

r 

axis  component,  p3,  of  the  warp  rate  vector,  p,  is  zero.  Then,  d3  can  also  be 
determined,  but  there  remains  the  indi stinguishabil ity  between  the  east-west 
component  of  earth  spin  vector  wander,  dl,  and  the  vertical  axis  component  of 
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crustal  warp  rotation,  a3 .  Bence,  inertial  neaaurements  don't  suffice  to 
define  the  deviation  of  the  earth  spin  vector,  nor  the  local  crustal  rotational 
warp.  Therefore,  we  cannot  use  inertial  measurements  alone  to  orient  to  the 
earth-fixed  reference  fraaie. 

RESOLUTION  OF  DISTINGUISHABILITJf 

Within  the  constraint  of  avoiding  the  use  of  astronomical  scale  fixes,  such  as 
lunar  laser  ranging,  star  shots,  etc.,  there  are  three  ways  of  dealing  with 
this  problem.  One  of  these  is  to  use  a  high-quality  theodolite,  affixed  to  the 
device  to  be  oriented,  and  landmark(s)  at  a  sufficient  distance  from  the 
viewing  point,  in  the  horizon  plane.  A  sufficient  distance  means  far  enough 
away  so  that  effects  of  landmark  crustal  rotation  are  attenuated  to  zero  (but 
not  so  far  as  to  give  excessive  refraction  errors)  .  How  far  remains  to  be 
determined,  as  a  function  of  what's  needed  and  what's  possible.  Clearly,  the 
wall  of  a  small  chamber  housing  the  apparatus  would  be  a  poor  landmark! 
distances  of  several  hundred  meters  would  probably  be  needed.  This  would 
enable  determination  of  crustal  rotation  and  rotation  rate  around  the  vertical 
axis  (a3  and  p3,  respectively),  to  the  accuracy  of  the  theodolite-operator 
sy  s  t  em . 

The  average  rotation  angle  around  the  vertical  would  be  the  average  of  some 
number  of  equi-time-spaced  readings  of  azimuthal  deviation,  taken  over  the 
observation  period.  The  average  rotation  rate  around  the  vertical  would,  in 
principal,  be  the  difference  between  an  initial  azimuthal  deviation  reading  and 
a  final  one,  divided  by  the  time  difference  between  the  two  readings.  In 
practice,  this  may  be  refined  by  best-fitting  to  all  the  azimuthal  deviation 
readings.  The  accuracy  of  individual  theodolite  fixes  is  about  1  arc-sec, 
rms.  If  this  is  random  independent  ''white  noise''  error,  100  successive  fixes 
would  average  the  error  down  to  0.1  arc-sec  rms.  Very  likely,  one  would  have 
to  take  successive  fixes  incessantly.  Well-correlated  fix  error  could  be  a 
major  problem. 


The  expected  plethora  of  theodolite  fixes  would,  of  course,  be  passed  into  a 
coaputer.  However,  this  method  could  still  be  very  labor-intensive,  because  of 
the  multitude  of  readings  to  be  taken,  and  because  of  the  indexing,  plunging, 
and  what-have-you  that  accompanies  theodolite  readings,  in  order  to  cancel  some 
of  the  systematic  errors.  Some  of  this  could  doubtless  be  automated,  but  there 
would  still  be  demands  on  the  operator.  For  these  reasons,  and  because  of  the 
eventual  environment  the  system  might  be  required  to  operate  in,  this  approach 
was  dropped  from  further  consideration. 

The  second  method  of  dealing  with  the  distinguishabil ity  problem  is  to  use 
independently  determined  earth  spin  vector  values.  Then,  we  know  the  values  of 
dl,  d2,  and  d3,  and  we  get  a3  and  p3  by: 

a3  =  (dl  -  el  +  R*gl*sin(L)  -  g *2) /(R*cos(L) ) 
p3  -  e3  +  R*g2*cos(L)  -d3 

The  earth  spin  vector  is  a  world  constant.  So,  it  can  be  determined  at  any 
convenient  location  and  its  value  then  transmitted  to  a  multitude  of  sites. 

The  value  is  then  easily  transformed  into  a  local  coordinate  system 
expression.  Specifically,  let  the  deviation  be  expressed  as  the  vectorial  sum 
of  a  polar  deviation  (i.e.,  a  deviation  in  magnitude) ,  DM,  a  deviation  along 
the  Greenwich  meridian,  DG,  and  a  deviation  along  the  90°  east  meridian,  DE. 
Then,  at  latitude  LN  north  and  longitude  LE  east,  the  locally-seen  earth  spin 
deviation  vector  is 

d  *  i*(DG*sin(LE)  +  DE*cos(LE)> 

+  j*(DM*cos(LN)  +  sin(LN) *(DG*cos(LE)  -  DE*sin(LE) ) ) 

+  k*(DM*sin(LN)  -  cos(LN) •<DG*cos(LE)  -  DE*sin(LE) ) ) 
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i.  e. , 


dl  -  DG*sin(LE)  +  DE*cos(LE) 

d2  -  DM*cos(LN)  +  sin(LN)*(DG*cos(LE)  -  DE*sin(LE) ) 

d3  »  DM*sin(LN)  -  cos(LN) •<DG*cos<LE)  -  DE*sin(LE) ) 

In  principle,  measuring  the  exact  earth  spin  vector  at  a  central  location  is  no 
easier  than  doing  it  at  a  multitude  of  sites.  But,  it  aay  be  possible  to  find 
locations  where  the  earth  crustal  warp  is  essentially  zero,  e.g.,  on  some 
special  underground  bedrock  granite  formation,  far  removed  from  seismic 
epicenters.  These  determinations  could  be  intermittently  checked  by  star 
sightings,  whenever  the  ''seeing' *  is  favorable.  The  assumption  is  that  some 
special  earth- spin-measurement  laboratory  can  be  set  up  at  some  favorable 
location(s)  to  provide  continuous  values  good  to  within  0.1  arc-sec  rms. 

There  is  a  third  possibility,  suggested  by  our  perception  of  the  ultimate 
application  by  the  customer  (AFGL) .  It  appears  that  the  system  could 
eventually  be  deployed  at  a  large  number  of  sites,  N,  with  measurements  taken 
fairly  frequently.  Furthermore,  it  would  be  reasonable  to  assume  that  the 
crustal  warps  at  all  N  sites  would  constitute  a  zero-mean  random  independent 
set.  Consequently,  the  multi-site  mean  crustal  rotation  around  vertical 
should  be  zero,  with  an  rms  deviation  around  that  value  of  magnitude  1/SQRT(N) 
times  the  rms  value  of  individual  site  crustal  rotations  around  the  multi-site 
mean.  If  so,  the  individual  determinations  of  site  crustal  rotation  plus  earth 
axis  deviation  should,  upon  multi-site  averaging,  produce  an  estimate  of  earth 
axis  deviation,  corrupted  by  1/SQRT(N)  times  the  rms  crustal  deviation. 
Actually,  the  computation  would  be  somewhat  more  complex  than  the  above, 
especially  if  the  sites  were  well  separated  on  the  earth  surface.  In  that 
case,  the  averaging  would  have  to  include  the  fact  that  earth  spin  axis 
deviation  will  be  seen  differently  at  each  of  the  N  locations.  It  would,  in 
fact,  be  a  straightforward  least-squares  best-fit  procedure. 
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DUAL-FREEDOM  STRAPDOWN  CONCEPT 


Northfinding/ leveling  is  performed  routinely  by  ell  inertial  navigation  systems 
as  a  preliminary  to  entering  the  navigation  phase.  Accuracy  of  northfinding  is 
typically  about  1  mrad,  achieved  in  15  min.  The  major  factors  affecting 
accuracy  are  gyro  biases  (fixed  and  random)  and  accelerometer  biases. 

In  order  to  enhance  performance,  ve  may  rotate  the  inertial  sensor  assembly 
(ISA)  about  one  or  more  axes,  called  carousel ing.  By  so  doing,  the  effect  of 
gyro  and  accelerometer  fixed  biases  may  be  partially  or  totally  eliminated, 
depending  on  the  number  of  axes  of  carousel ing. 

In  a  pure  strapdown  system  (see  Figure  la),  attitude  error  is  essentially  the 
time  integral  of  all  gyro  error  terms.  In  order  to  reduce  the  error,  the 

inertial  sensors  may  be  mounted  to  a  rotating  base,  rather  than  to  a  fixed  one 

(see  Figure  lb).  Partial  commutation  or  averaging  of  systematic  gyro  errors 
results.  There  is  a  reduced  attitude  error  for  axes  in  the  plane  of  rotation, 

but  none  for  the  the  out-of-plane  one.  (It  does  not  help  to  dispose  the  three 

gyros  so  that  all  three  gyro  axes  partake  of  the  rotation.) 

Further  attenuation  is  provided  by  a  dual-freedom  strapdown  approach,  in  which 
the  plane  of  the  rotating  turntable  is  itself  being  rotated  (see  Figure  lc). 

If  both  rotational  axes  are  properly  controlled,  it  is  possible  to  get  almost 
bounded  attitude  errors  over  very  long  time  periods,  due  to  constant  gyro 
errors  such  as  bias,  scale  factor  error,  and  misalignment.  There  is,  however, 
no  reduction  of  the  error  due  to  dithered  laser  gyro  random  walk. 

Having  thus  dealt  with  systematic  error  sources,  the  final  step  is  to  rate-bias 
the  gyros  to  eliminate  random  walk.  The  three  gyros  are  maintained  in 
uninterrupted  unidirectional  rotation,  rather  than  being  dithered  in  periodic 
motion.  This  eliminates  the  lock-in  error  increments  that  occur  as  the  gyro 
rotation  rate  goes  through  zero  twice  per  dither  cycle,  and  hence  eliminates 
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Figure  1.  Evolution  of  Dual-Freedom  Strapdown  Triad 


the  resultant  ran  don  walk.  To  do  this  on  the  daal-freedoai  approach,  we  elevate 
the  gyro  input  axes  above  the  equatorial  plane,  so  that  all  three  gyros  receive 
a  portion  of  the  angular  rate  about  the  inner  axis,  and  we  slip-ring  the  inner 
axis,  to  allow  continuous  rotation. 

To  suaunarize,  pure  strapdown  is  the  simplest  approach,  but  has  the  greatest 
attitude  error  growth.  The  single-  and  dual-freedom  approaches  attenuate, 
respectively,  some  and  all  of  the  attitude  error  due  to  fixed  gyro  errors, 
without  any  reduction  of  that  part  due  to  random  walk.  Finally,  dual-freedom 
strapdown  with  rate-biasing  eliminates  virtually  all  errors.  This  is  the 
configuration  that  we  have  selected  for  analysis  and  simulation. 

DUAL-FREEDOM  STRAPDOWN  ROTATIONAL  STRATEGY 

Appendix  A  shows  that  two  conditions  suffice  to  enable  successful  dual-freedom 
strapdown  rotational  strategy 

1.  The  inner  axis  rate  is  an  integer  multiple  of  the  outer  axis  rate. 

2.  The  outer  axis  reverses  direction  at  360°  traverse  intervals. 

These  ensure  the  commutation  of  gyro  and  accelerometer  bias,  scale  factor,  and 
misalignment  errors. 

The  Appendix  A  analysis  assumed  that  the  outer  axis  can  reverse  direction 
instantaneously,  and  that  there  are  no  nonlinearity  errors  in  the  sensor 
outputs.  But,  the  analysis  can  be  extended  by  expressing  all  quantities 
appearing  in  Appendix  A  tables  as  Fourier  series  over  one  cycle  of  the  outer 
axis.  This  leads  to  a  refined  dual-axis  rotational  profile  that  gives  periodic 
cancellation  of  the  following  sensor  errors. 


*  gyro  bias  errors 

*  gyro  linear  scale  factor  errors 

*  gyro  misalignment  errors 

*  gyro  scale  factor  asymmetry  and  nonlinearity  errors 

*  accelerometer  bias  errors 

*  accelerometer  linear  scale  factor  errors 

*  some  accelerometer  scale-factor  nonlinearity  errors  (all  odd  powers,  Oth 
and  2nd  even  powers  only) 

*  accelerometer  misalignment  errors 

An  example  of  error-cancelling  gimbal  angle  profiles  is  shown  in  Figure  2.  The 
outer  gimbal  rates  consist  of  zero  rate  segments,  constant  rate  segments,  and 
ramp  rate  segments.  More  generally,  error  cancellation  is  provided  when  each 
segment  of  the  outer  axis  rate  profile  is  symmetric  about  its  midpoint. 


Figure  2.  Dual-Freedom  Sensor  Assembly  Rotational  Profile 


SECTION  3 


SENSOR  AND  SYSTEM  ERROR  MODELS 

GYRO  RATE  WHITE  NOISE 

The  effect  of  rate  white  noise  on  system  attitude  error  can  be  defined  by  the 
simpl e  model : 

Cgng 
■  CAng 

where  C  and  A  are  defined  by: 

C  =  transformation  matrix  from  sensor  reference  axes  to  local-vertical 
reference  axes 

A  =  transformation  matrix  from  gyro  axes  to  sensor  reference  axes 

The  covariance  matrix  of  V.  for  equal  intensity  white  noise  in  each  gyro  is 
governed  by: 

p  =  caaTCtn 

where: 

P  =  E('P'irT) 

and  N  is  the  gyro  white  noise  spectral  intensity. 
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If,  as  is  frequently  the  case,  the  matrix  A  is  orthogonal,  and  since  C  is 
orthogonal  by  definition,  the  covariance  equation  reduces  to: 

P  =  N 

or: 

P  =  Nt, 

which  is  to  say: 


where  t  is  the  elapsed  time.  Therefore,  rotation  of  the  sensor  assembly 
(implicit  in  C)  will  do  nothing  to  reduce  the  effect  of  the  gyro  white  noise. 

SENSOR  QUANTIZATION 

A  thorough  assessment  of  sensor  quantization  errors  is  given  in  Appendix  C, 
which  defines  a  covariance-based  model  suitable  for  performance  analysis 
purposes.  The  essence  of  the  quantization-induced  error  can  be  described  more 
simply  by  considering  a  single-axis  rotational  system.  One  finds  that  the  gyro 
quantization  error  leads  to  an  attitude  random  walk  error,  which  is  defined  by: 

Oq,  =  uoq(T/2)  l/2tl/2  (3-2) 

where: 

<jy  *  standard  deviation  of  attitude  error 

to  *  sensor  assembly  rotational  rate 

Oq  *  standard  deviation  of  the  gyro  quantization  error 


T  “  attitude  update  interval 
t  *  time 

and: 

®q  -  P/'JlT 

where  P  is  the  gyro  pulse  size. 

A  simplified  analysis  of  the  effect  of  accelerometer  quantization  leads  to  a 
formula  analogous  to  Equation  3-2  that  defines  a  random  walk  error  growth  in 
velocity. 

It  is  clear  from  Equation  3-2  that  the  random  walk  error  is  directly 
proportional  to  rotation  rate,  and  hence,  is  a  major  error  contributor  for 
relatively  high  rotation  rates. 

PARAMETER  VARIATION 

The  error  cancellation  properties  associated  with  any  rotational  strategy  are 
applicable  only  when  the  sensor  errors  are  constant.  In  reality  the  sensor 
biases,  misalignments,  etc.,  will  randomly  vary  due  to  environmental  and  other 
causes . 

Three  types  of  random  sensor  error  variations  are  considered  below. 

1.  First-order  Gauss-Markov  process,  defined  by  the  first-order  stochastic 
differential  equation: 

k  +  Pjk  =  q 
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where: 


k  ■  random  variation  in  the  sensor  error  coefficient 
n  ■  white  noise 

01  *  1/correlation  time  of  Gauss-Markov  process 

The  power  spectral  density  of  the  white  noise  input  is  chosen  according 
to: 


%  -  2p!02 

where : 

Nj  91  power  spectral  density  of  the  white  noise  input 
a  -  standard  deviation  of  the  steady-state  sensor  error  variation 

A  typical  sample  of  a  first-order  Gauss-Markov  time  history  is  shown  in 
Figure  3a.  The  important  feature  is  that,  even  for  a  one-hour 
correlation  time,  the  response  shows  significant  high  frequency 
activity. 

2.  Second-order  Gauss-Markov  process,  defined  by: 
k  +  202k  +  0*k  "  1 

where  k  is  the  error  parameter  variation  and 
02  =  2.146/TC 


v.  • 


.-V 
■  -o  c 


v"  It'  » 
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rm  (ml 

a.  First-Ordar  Ganss-Markov  (o  *  1  unit,  t  » 


THM  (ml 


b.  Second-Order  Gaoss-Narkov  (or  ■  1  unit,  t  = 


mac  (mi 

c.  Random  Walk  (o^*  “  1  unit) 


1  hr) 


1  hr) 


Fignre  3.  Sample  Time  Histories  of  Three  Stochastic  Processes 


in  which 


Tc  ■  correlation  ti»e  of  tha  aacond-ordar  Ganas-Markov  process 

The  apactral  danaity  of  tha  white  noiaa  function,  q,  ia  defined  by 
N2  -  4032oJ 

where,  as  before,  or  ia  the  atandard  deviation  of  the  ateady-state  sensor 
error  variation. 

A  typical  tine  history  for  a  second-order  Ganas-Markov  proce:  is  shown 
in  Figure  3b.  Note  that,  with  one-hoar  correlation  tine,  the  response 
is  essentially  free  of  all  of  the  high  frequency  variations  that  are  so 
prominent  in  the  first-order  Ganas-Markov  proceas  with  the  same 
correlation  time. 

3.  A  third  type  of  random  sensor  error  parameter  variation  is  defined  by: 

k  =  q 

which  constitutes  a  random  walk  variation  in  the  parameter.  A  sample 
time  history  for  a  random  walk  parameter  variation  is  shown  in 
Figure  3c.  The  spectral  intensity  of  the  white  noise  function,  q,  is 
defined  by: 

N  ■  <®hr/60>2 

where  is  the  standard  deviation  of  the  parameter  growth  at  one  hour. 
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COMPLIANCE-RELATED  ERRORS 


In  a  dual-freedom  strapdown  system,  rectified  attitude  drift  may  generated  by 
g-induced  gyro  input  axis  compliance  effects. 

Attitude  error  drift  rate  due  to  gyro  input  axis  compliance  may  be  characterized 
by: 

S  "  C8N4  (3-3) 


where : 

*  attitude  error  vector  (with  components  along  nor th- east-down  (NED)) 

C  =  transformation  matrix  from  a  set  of  orthogonal  sensor  reference  axes  to 
the  NED  frame 

6M  =  error  in  the  transformation  matrix  from  gyro  input  axes  to  the  sensor 
reference  frame 

w*  =  nominal  angular  rate  vector  measured  by  the  gyro  triad 


The  nominal  angular  rate  vector  measured  by  the  gyros  is  related  to  the  nominal 
angular  rate  vector  coordinatized  in  the  orthogonal  sensor  reference  axes 
according  to: 


8  U~1  r 

(1)  =  M  (i) 

~0  O  - o 

where : 


(3-4) 


u»r  =  nominal  angular  rate  vector  coordinatized  in  orthogonal  sensor 
o 

reference  axes 


M  *  =  nominal  transformation  matrix  from  orthogonal  sensor 
o 

reference  axes  to  gyro  input  axes 
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therefore 


6 MM  1  +  M  6M_1M  =  0 

o  o  o 

and 

6M  =  -M  8M-1M  (3-5) 

o  o  v 

Substituting  Equation  3-4  and  Equation  3-5  in  Equation  3-3  gives 
•  - 

V  *  -CM  6M  «r  (3-6) 

—  o  ~o 

which  is  the  basic  attitude  error  equation. 

To  complete  Equation  3-6,  the  following  equations  defining  C  Mad  gr  are  defined. 


CaCp 

-CaSp 

-• 

Sa 

c  = 

sp 

cp 

0 

(3-7) 

-SaCp 

SaSp 

Ca 

r  T 

uo  =  (asinpt  acosflt  (1)  (3-8) 

where : 

a  s  outer  gimbal  angle 
8  =  inner  gimbal  angle 


For  a  gyro  triad  whose  axes  are  nominally  orthogonal  and  coincident  with  the 
sensor  assembly  reference  axes: 


Mo  *  I 

and  the  error  in  M0-1  can  be  expressed  in  general  by: 


(3-9) 


8M~1  = 


0 

8aj 

6b! 

8  82 

0 

6b2 

8  *3 

8b3 

0 

(3-10) 


where  6a^  and  6bj  are  the  gyro  inpnt-axis  misalignment  angles. 


Now  assume  that  each  misalignment  angle  error  can  be  defined  in  terms  of  the  thre 
components  of  gravity  in  the  orthogonal  sensor  reference  frame  according  to: 


.  i  .  .r 
8a  *  k 


11*1  +  ki2*2  +  ki3“3 


r  +  k a5 


8bi  =  milAl  +  mi2A7  + 


i3  3 


(3-11) 

(3-12) 


where  the  kjj  and  m^j  are  compliance  coefficients,  and  A*,  A*,  and  A|  are  the 
components  of  the  gravity  reaction  acceleration,  defined  from  Equation  3-7  by: 


Ar  -  -g(-SaCp  SaSp  Ca)T 
where  g  is  the  acceleration  dne  to  gravity. 


Substituting  Equations  3-7,  3-8,  3-9,  3-10,  3-11,  and  3-12  into  Equation  3-6 
yields  an  attitude  error  equation  that  takes  the  following  form. 


?  *  f  (flj  P#  A*  Pi  kj j »  kj 2  •  •  •  *  1 »  2*  •  •  *  * 

where: 

?  -  Wn  'J'e  ^d)1 

•  ® 

and  all  terns  in  the  linear  expansion  have  either  a  or  0  as  a  multiplier. 


Each  component  of  $  can  be  determined  at  the  end  of  a  complete  gimbal  period  by 
direct  integration  of  In  the  process  of  doing  this,  one  can  easily  eliminate 
terms  involving  a  by  inspection.  For  example,  a  term  such  as 

I  I  ® 

J  asinatcosfjt  dt 
o 


is  identically  zero  at  the  completion  of  every  outer  gimbal  360°  traversal  (of 
duration  T) ,  while  a  term  such  as 

/* P  sin2 p  dt 


has  a  non-zero  mean  value  at  the  end  of  each  outer  axis  period. 


SECTION  4 


SYSTEM  ERBOR  STUDIES 

ERROR  PROPAGATION  IN  DUAL-FREEDOM  SYSTEM 

In  a  rotating  system,  error  propagation  depends  on  the  nature  of  the  random 
sensor  errors.  For  example,  consider  the  simplified  block  diagram  shown  in 
Figure  4  for  attitude  error  propagation  in  a  single-freedom  system.  The 
variance  of  the  attitude  error  can  be  determined  from: 

o2(t)  =  N /*w  ^  t,x) dt 
l||  o 

where: 

w(t,x)  =  response  of  the  system  at  time  t  to  an  impulse  applied  at  time  x 

N  =  power  spectral  density  of  the  white  noise  input 

<j2  =  variance  of  the  attitude  error 

V 

For  the  first-order  Gauss-Markov  process,  the  impulse  reponse  function  is  given 
by: 

w(t,x)  =  /le  ^1^*  T^sin«i)xdx 
o 

and  for  the  second-order  Gauss-Markov  process: 

w(t,x)  *  /t(t-x)e  ^2**  T*sinwxdx 
o 

For  the  random  walk  process: 

w(t,x)  *  /*  U(x)  sin  «xdx 
o 

where  U  denotes  the  unit  step  function. 
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a.  First-order  Gauss-Markov  process 


b.  Second-order  Ganss-Markov  process 


c.  Random-walk  process 


T)  =  UNIT  INTENSITY  WHITE  NOISE 

a  =>  STANDARD  DEVIATION  OF  GYRO  BIAS  VARIATIONS  IN  STEADY  STATE  (GAUSS-MARKOV  PROCESS) 

gj  =  ROTATION  FREQUENCY 

i j)  =  ATTITUDE  ERROR 

P ,  =  (/(CORRELATION  TIME  OF  FIRST-ORDER  GAUSS-MARKOV  PROCESS) 

P  2  =  2. 146/(CORR ELATION  TIME  OF  SECOND-OROER  GAUSS-MARKOV  PROCESS) 

ohr  =  STANDARD  DEVIATION  OF  GYRO  BIAS  VARIATION  AT  1  hr  (RANDOM  WALK  PROCESS) 

Figure  4.  Attitude  Error  Propagation  in  a  Rotational  System 


The  asymptotic  characteristics  of  the  solutions  can  be  established  by  carrying 
out  the  integrations  analytically  and  discarding  all  terms  in  the  solution  that 
decay  with  time.  The  result  is  that  the  attitude  error  variances  become 
asymptotically  equal  to: 

a*  »  0^(u»)t/2  (first-order  Gauss-Markov)  (4-1) 

2 

a  =  0^  (<«»)  t/2  (second-order  Gauss-Markov)  (4-2) 


where  <>1<  w)  is  the  power  spectral  density  of  the  first-order  Gauss-Markov 
process  at  the  rotational  frequency,  w,  and  is  defined  by: 


and  02(<u)  is  the  power  spectral  density  of  the  second-order  Gauss-Markov 
process  at  the  rotational  frequency,  and  is  given  by: 


02U) 


.tt3  2 
*2m 


<P*  ♦  »2)2 


For  the  random  walk  parameter  variation 


a  =  <M«)t 


where: 


®r(“’  ’  <V/60|2/“2 


(4-3 


It  is  seen,  therefore,  that  all  three  processes  lead  to  a  random  walk  error 
propagation  tendency  in  attitude.  It  is  also  clear  that,  for  the  Gauss-Markov 
processes,  the  variance  of  the  random-walk  component  of  attitude  error  is 
directly  proportional  to  the  amount  of  power  in  the  process  at  the  rotational 
frequency,  <■>.  Comparing  and  $2‘  the  ratio  of  powers  at  the  rotational 
frequency  is  determined  as 


fC  +  <* 


.  2.2 
+  u  ) 


p^p2**.2)2 

2p2(p2-H„2) 


For  example,  consider: 


Px  *  1/1800 
P2  *  2.146/1800 
w  =  2n/60 


which  yields  the  resol t  that 
02  (<a)  (ec)  =  1799 

Therefore,  the  standard  deviation  of  the  random  walk  component  is  greater  by  a 
factor  of  ^1799  =  42.4  for  the  first-order  Ganss-Markov  process  compared  to  the 
second-order  Gauss-Markov  process  for  the  example  chosen. 

It  is  clear  from  Equations  4-1,  4-2,  and  4-3  that  each  type  of  parameter 
variation  leads  to  a  random-walk  in  attitude.  It  can  be  shown  that,  for  a 
Gauss-Markov  process  of  any  order  whatsoever,  the  attitude  error  growth  is 
defined  by  the  general  form: 

<4  =  f>N(<i>)t/2 

where  0jv|(<t>)  is  the  power  spectral  density  of  the  nt*1-order  Gauss-Markov  process 
at  the  rotational  frequency  ui.  The  implication  is  that  the  choice  of  a  random 
process  for  representing  sensor  variations  is  not  critical.  It  is  important 
only  that  the  power  at  the  rotational  frequency  have  some  meaningful  value. 

This  is  true  for  performance  analyses  and  for  Kalman  filter  design  and 
evaluation. 

COVARIANCE  ANALYSIS  AND  KALMAN  ESTIMATION 

Consider  a  dynamic  system  whose  behavior  is  defined  by  the  following  set  of 
discrete  linear  equations: 

=  Pn  ^n-1  +  ^n’ln  (4—4) 

where : 

X  =  vector  of  states 

q  =  vector  of  random  (zero  mean)  noise  sequences 
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0n  =  state  transition  matrix  from  (n-l)tl1  to  n^  iterations 
Bn  =  noise  distribtntion  matrix  at  the  n iteration 


For  a  given  0  and  B,  the  state  X  will  have  a  time  variation  determined  by  the 
particular  noise  sequence,  q,  and  the  initial  condition,  X0,  which  in  general 
is  taken  to  be  a  randomly  distributed  quantity.  Since  the  noise  sequence,  q, 
has  an  infinite  number  of  realizations  and  the  initial  condition  error  can 
assume  an  infinite  number  of  values,  the  system  of  equations  given  by  Equation 
4-4  has  an  infinite  number  of  solutions.  Because  this  is  true,  attention  must 
be  focused  on  the  statistical  behavior  of  Equation  4-4  rather  than  on  specific 
solutions. 


The  most  natural  and  useful  way  of  characterizing  the  behavior  of  Equation  4-4 
is  to  compute  the  statistical  parameters  that  define  the  bounds  on  the  state 
vector,  X.  The  statistical  bounds  on  the  components  of  X  are  found  by  solving 
the  covariance  matrix  equation  associated  with  Equation  4-4,  which  takes  the 
recursive  form: 

pn  “  ®n  Pn-1  ®nT  +  Bn  <*n  BnX  (4-5 

where  P  is  the  state  covariance  matrix  of  the  state  vector,  X,  defined 
explicitly  by: 

P  *  tPii] 


Pij  E(XjXj) 


3 


in  which  E  denotes  the  expectation  operator.  It  is  seen  that  the  individual 
variances  of  the  components  of  X  are  defined  by  the  diagonal  elements  of  P, 
with  the  joint  expectancies  (cross-correlations)  being  defined  by  the 
off-diagonal  elements  of  P.  The  matrix  Q  in  Equation  4-5  is  the  covariance 
matrix  of  the  driving  noise  vector,  q,  which  is  defined  by: 

Q  =  Uijl 

and 

qij  =  E(niT|j) 

Suppose  that  the  discrete  process  defined  by  Equation  4-4  represents  the  true 
dynamic  propagation  characteristics  of  a  given  linear  system.  Suppose  further 
that  at  some  time  a  measurement  is  made  (using  an  external  measuring  device) 
that  allows  a  specific  linear  combination  of  the  states  to  be  directly 
monitored.  A  generalized  way  of  stating  this  in  mathematical  form  is 

Ym  *  HmX  +  (4-6) 

where : 

Ym  =  vector  of  measurements  at  time  t,„ 

Hm  =  measurement  matrix  at  time  tg, 

=  measurement  noise  vector  applicable  to  m^E  measurement 

and  it  is  assumed,  in  the  general  case,  that  a  number  of  independent 
measurements  may  become  available  simultaneously. 

Consider  now  the  question  of  how  one  could  use  the  information  introduced 
through  a  series  of  measurements  of  the  form  given  by  Equation  4-6  to  optimally 
estimate  the  state  X  in  a  sequential  fashion.  This  is  the  central  problem 
addressed  by  Kalman  estimation  theory,  and  has  the  following  solution. 


After  each  measurement  (of  a  sequence  of  measurements)  is  made,  the  estimate  of 
the  state  X  is  refreshed  by  the  two-step  procedure: 


Ut  )  =  0(t  ,t  At  >  (4-7) 

m  m  m-1  m- 1 

$(t  )  =  $(t  >  +  K  (Y  -  H jT(t  ))  (4-8) 

n  in  m  in  n  id 


where: 


th 


X  (t  )  =  best  estimate  of  vector  X  just  prior  to  the  m  measurement 

ID 


$(t  )  =  best  estimate  of  vector  X  after  incorporating  the  information 


th 


introduced  by  the  m  measurement 


th 


0(t  ,t  )  -  state  transition  matrix  spanning  the  time  between  the  m  and 
m  m-1 

(m-1)  measurements 


Kg,  =  Kalman  gain  matrix  at  the  m(^  measurement  point 
Ym  =  vector  of  measurements 

Hm  =  measurement  (observation)  matrix  at  m**1  measurement 


with  Kg,  defined  by: 


K  =PHT(HPHT+R)1 
id  mm  id  m  m  nt 


(4-9) 


in  which 

P  -  a  priori  error  covariance  matrix  of  the  vector  X 
m 

T 

R  3  measurement  error  covariance  matrix  =  E($  £  ) 
m  mm 

and  the  a  priori  error  covariance  matrix  Pq  is  computed  by  using  Equation  4-5 
over  the  period  from  the  to  (m-1) 111  measurements. 


After  processing  of  the  m**1  measurement,  the  error  covariance  matrix  of  the 
state  X  is  modified  to  reflect  the  benefit  of  incorporating  the  nev  information 
introduced  by  the  measurement,  as  follows. 

P  -  (I  -  K  H  )P_  (4-10) 

m  m  m  m 

and  Pm  is  referred  to  as  the  4a  posteriori"  error  covariance  matrix.  The  form 
given  by  Equation  4-10  is  applicable  when  the  Kalman  filter  is  fully  optimal* 
that  is,  the  filter  is  a  full-state  filter  in  which  all  components  of  X  are 
fully  acounted  for  in  the  mathematic  model,  and  further,  are  re-estimated  after 
each  successive  measurement  is  made  available. 


Formulation  of  Dual-Freedom  System  Simulation 


The  northfinding  process  is  assumed  to  take  place  in  real-time  (rather  than  by 
data  post  processing)  via  optimal  estimation  (Kalman  filtering)  techniques. 

The  alignment  Ealman  filter  is  based  on  the  fact  that  the  reference  origin  of 
the  sensor  assembly  has  zero  velocity.  Then,  the  velocity  indicated  by  the 
system  becomes  a  measure  of  the  alignment  errors  (azimuth  and  level). 


A  discrete  error  equation  for  expressing  deviations  in  attitude  and  velocity 
for  the  northfinding  mode  is  given  as  (see  Appendix  B) : 


V  *11  *n-l  +  *12  ~n-l 


(4-11) 


where  Y  is  the  vector  of  attitude  and  velocity  error  states  defined  by: 


Y  =  (V  &V)T 


(4-12) 


in  which : 


V  =  attitude  error  vector 
6V  =  velocity  error  vector 
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and  X  is  the  vector  of  sensor  errors,  which  are  represented  in  general  by 
second-order  Gauss-Markov  processes: 

X  =  (Xj  X2  ...  (4-13) 

Each  component  of  the  sensor  bias  vector,  is  defined  by  two  first-order 
differential  equations  according  to  (see  Figure  4b): 

Xi  +  P2*i  =  Ui  (4"14) 

Ui  +  P2Ui  *J4p2^n  (4'15) 

where  q  is  a  unit  intensity  white  noise  function,  and 
02  ~  2 . 146/ (correlation  time) 

The  vector  X  represents  the  short-term  instability  in  the  sensor  error  parameter 
vector.  Only  the  short-term  instability  is  of  importance.  The  constant 
components  of  sensor  error  are  dealt  with  by  the  error  commutation  properties 
inherent  in  the  sensor  assembly  rotational  profile  (see  Figure  2).  The  single 
exception  is  vertical  drift  rate,  due  to  input-axis  compliance.  Since  this 
effect  is  highly  stable  over  time  (based  on  experimental  evidence),  it  is 
assumed  that,  prior  to  northfinding,  the  vertical  drift  rate  has  been  calibrated 
and  compensation  provided.  Alternatively,  the  vertical  drift  rate  may  be 
fine-calibrated  during  the  northfinding  process  by  adding  an  additional  state  to 
the  Kalman  filter.  In  either  case,  the  effect  of  the  vertical  drift  rate  can  be 
made  negligibly  small. 

Equations  4-11  through  4-15  describe  the  evolution  of  the  error  state.  The 
discrete  uncorrelated  process  noise  vector,  £,  is  defined  to  include  gyro  random 
walk  and  sensor  quantization  errors. 


Hie  measurement  equations  for  the  northfinding  prooess  are  defined  by: 


yi  -  SVi 

♦  »*1 

(4-16) 

72  -  &V2 

+  1*2 

(4-17) 

where: 

y^,  “  indicated  velocity  ontpnts  of  the  navigation  system  along  north 

and  east  axes,  respectively 

6V^,  6 =  errors  in  the  north  and  east  velocity  components,  respectively 

|i^,  =  measurement  errors  (spurions  velocity  components  due  to  random 

base  motion  effects) 

Based  on  Equations  4-11  through  4-17,  an  alignment  Kalman  filter  is  defined.  In 
the  present  study,  the  Kalman  filter  update  interval  is  chosen  to  be  the  same  as 
the  period  of  the  outer  axis  rotational  profile  (4  min) .  By  so  doing,  the 
achievable  alignment  accuracy  becomes  independent  of  any  latent  sensor 
calibration  errors.  Only  short-term  instabilities  in  the  sensor  errors  (as 
defined  by  Equations  4-14  and  4-15)  affect  the  accuracy  of  the  alignment.  By 
enforcing  this  constraint,  the  computational  requirements  of  the  alignment 
Kalman  filter  are  greatly  reduced,  without  loss  of  performance. 

A  set  of  la  sensor  errors  that  was  assumed  for  assessing  the  performance  of  the 
three-axis  navigator  approach  is  defined  in  Table  1.  The  gyro  random  walk 
coefficient  of  .45  x  10“ 4  deg  at  1  hr  is  based  on  experimental  results  obtained 
for  the  Honeywell  GG1389  laser  gyro.  The  quantization  lo  errors  are  based  on 
the  assumption  of  a  uniform  distribution  for  the  quantization  error,  and  the  use 
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TABLE  1 .  la  ERRORS  ASSUMED  FOR  THREE- AX IS  SIMULATION 


Random  walk  =  0.00004 5 °/rt-hr 

Quantization  =  0.022  arc-sec 

Gyro  bias  short-term  variation  =  0.01°/hr 

Gyro  scale  factor  short-term  variation  =  0.05  ppm 

Gyro  inpnt  axis  misalignment  short-term  variation  =  0.05 


arc-sec 


Accelerometer 


Quantization  =  0.0006  ft/sec 
Accel  bias  short-term  variation  >*  0.10  micro-g 
Accel  scale  factor  short-term  variation  =  0.5  ppm 
Accel  input  axis  misalignment  short-term  variation 


arc-sec 


Other 

Short-term  translational  base  motion  -  0.005  ft/sec 


of  a  phase-lock  loop  to  increase  the  resolution  of  the  gyro  output  to  0.075 


arc-sec . 


The  short-term  sensor  errors  are  assumed  to  be  second-order  Gauss-Markov 
processes  with  a  correlation  time  of  1  hr.  However,  the  lo  values  given  in 
Table  1  are  not  meant  to  imply  bounded  behavior.  Rather,  the  intent  is  to 
characterize  their  short-term  variation  over  times  comparable  to  the  period  of 
the  outer  axis  (4  min  in  this  case). 


The  three-axis  system  covariance  and  error  sensitivity  analyses  were  first 
framed  around  the  assumption  of  using  three  GG1300  RLGs,  with  a  complete  set  of 
rms  error  parameter  values  assumed  as  below.  Errors  are  assumed  random 
independent,  and  with  Gaussian  distributions,  except  for  quantization  errors, 
which  are  assumed  to  be  sum-bounded,  and  to  have  triangular  distributions. 

Gyro:  Honeywell  GG1300  RLGs,  in  rate-biased  mode  (i.e.,  prolonged 

unidirectional  rotation) 

Angle  random  walk 
Quantization/spillover 
Bias  random  walk 
Scale  factor  random  walk 

Accelerometers:  3  Bell  11s,  or  equivalent 

Pulse-weight  quantization  0.0006  ft/sec 

Acceleration  random  walk  0.1  micro-g’s/rt-hr 

Scale  factor  random  walk  0.5  ppm/rt-hr 

Input  axis  instability  0.5  arc-sec/rt-hr 

Bias  errors  are  cancelled  by  carouseling  around  two  indexation  axes,  over  a 
4-min  cycle.  Consequently,  the  system  puts  out  updates  of  three-axis  attitude 
15  times  per  hour.  At  each  update,  the  system  output  is  reset,  by  ’’reminding" 
the  system  computations  that  the  system  has  undergone  zero  displacement  over  the 
preceding  4  min,  and  has  zero  velocity. 


0.0001°/rt-hr 
0.03  arc-sec 
0.01°/hr-rt-hr 
0.05  ppm/rt-hr 


The  reset  process  is  clearly  seen  in  Figure  5.  which  shows  the  evolution  of  rms 
drift  error  around  the  north-south  axis,  as  produced  by  computer  simulations. 
The  error  locus  is  the  lower  envelope,  with  an  asymptotic  value  of  about  0.085 
arc-sec.  The  simulation  is  quite  similar  for  east-west  axis  drift.  In 
Figure  6,  we  see  the  simulation  of  vertical  axis  drift  error.  This  has  an 
apparent  asymptote  of  about  1.4  arc-sec.  The  drift  is  smaller  around  the 
horizontal  axes  because  the  gyros  are  greatly  helped  by  the  tiltmeter  aspect  of 
the  accels.  Horizontal  axis  drift  is  dominated  by  the  accel  quality  and 
vertical  axis  drift  by  the  gyro  quality. 

The  Honeywell  GG1300  RLG,  on  which  the  foregoing  results  were  obtained,  used  to 
be  our  best  gyro.  Now,  we  have  the  highrperf ormance  GG1389  RLG,  with  the 
following  characteristics: 


SIMULATION  OF  NORTH-SOUTH  rms  ERROR 
OF  3-AXIS  NAVIGATOR,  USING  GG1300s, 
SAMPLED  AT  4-min  INTERVALS 


TIME  Otrs) 

Figure  5.  Simulated  North-South  rms  Error,  Using  GG1300s 


SIMULATION  OF  VERTICAL  AXIS  rim  ERROR 
OF  3-AXIS  NAVIGATOR,  USING  GG1300t, 
SAMPLED  AT  4-min  INTERVALS 
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Figure  6.  Simulated  Vertical  zma  Error,  Using  GG1300s 


Gyro:  Honeywell  GG1389  RLG.  in  rate-biased  mode  (i.e.  prolonged 

unindirectional  rotation) 


Angle  random  walk 
Quantization/ spill  over 
Bias  random  walk 
Scale  factor  random  walk 
Input  axis  instability 


0. 000045°/ rt-br 
0.02  arc-sec 
0.01°/hr-rt-hr 
0.05  ppm/rt-hr 
0.05  arc-sec/ rt-hr 


Ve  see  from  Figures  7  and  8  that  the  simulations  show  drift  errors  around  the 
east-west  and  north-south  axes  of  about  0.075  arc-sec.  This  is  about  the  same 
as  with  GG1300s.  As  previously  mentioned,  this  part  of  the  drift  is 
accel-limited.  In  Figure  9  we  see  that  drift  around  the  vertical  is 
asymptotically  0.55  arc-sec;  as  could  be  expected,  this  gyro-limited  term  is 
improved  by  going  from  the  GG1300  to  the  GG1389.  However,  it  is  still  too  high 
for  a  subarc-second  system. 


Discussion  of  Simalation  Results 

The  reason  that  the  vertical  axis  drift  is  asymptotically  bounded  (rather  than 
being  an  unbounded  random  drift  to  match  the  unbounded  random  walk  of  the 
vertical  component  of  the  gyro  outputs),  is  as  follows:  The  vertical  axis  drift 
attempts  to  ''tell''  the  system  that  the  horizontal  gyro  sensing  plane  is  in 
unbounded  rotation.  If  this  were  really  so,  then  the  direction  of  the 
horizontal  component  of  earth  rate,  as  seen  by  the  system,  should  also  appear  to 
be  unboundedly  rotated.  When  this  fails  to  occur,  the  system  starts  to 
''ignore''  the  vertical  axis  drift.  The  final  result  is  a  compromise  between 
vertical  axis  random  angle  drift  and  horizontal  plane  angular  rate  white  noise. 


0.«  p 

SIMULATION  OF  EAST-WEST  rim  ERROR  OF 
3-AXIS  NAVIGATOR.  USING  GG1389s, 
SAMPLED  AT  4  MIN  INTERVALS 


Figure  8.  Simulated  North-South  ms  Error,  Csing  GG1389s 
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Figure  9.  Simulated  Vertical  ms  Error,  Using  GG1389s 


The  azimuth  error  is  estimated  via  the  two  horizontal  velocity  components,  and 
this  process  tends  to  follow  the  law 


N1'*  .-1/2 

u>  cX 


where 


(4-18) 


N  =  intensity  of  gyro  white  noise  drift 

<de  =  earth's  rotational  rate  (0.729  x  10“*  rad/sec) 

cX  -  cosine  of  latitude 

=  lcr  error  in  estimated  azimuth 

Simultaneously,  the  azimuth  error  tends  to  grow  due  to  the  component  of  total 
gyro  random  drift  acting  in  the  vertical  direction.  This  increase  in  azimuth 
error  follows  the  law 


Nl/2  tl/2 


(4-19) 


Therefore,  the  decay  of  error,  as  predicted  by  Equation  4-18,  is  equal  to  the 
growth  of  error  predicted  by  Equation  4-19  when: 


1 


w  cX  t 
e 


1/2 


t 


1/2 


or  when 

t  * 


w  cX 
e 


For  X  *  45°,  this  gives: 

t  *  19402  sec  =  5.4  hr 
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This  is  reasonably  close  to  the  azimuth  error  convergence  time  constant  seen  in 
Figure  9.  The  basic  reason  we  see  this  irreducible  azimuth  error  is  that  the 
overall  model  included  the  possibility  of  azimuth  drift  over  the  measurement 
period. 

If  one  assumes  that  the  system  maintains  an  essentially  fixed  orientation 
relative  to  north  over  an  indefinite  time,  further  accuracy  of  northfinding  is 
possible.  A  simple  approach  would  be  to  average  a  series  of  northfinding 
results,  each  obtained  over  some  fairly  short  period  of  time  (say  about 
1/4  hr).  Then,  the  composite  accuracy  would  be  improved  by  a  factor  equal  to 
the  square  root  of  the  total  measurement  time.  If,  as  an  example,  64 
northfinding  experiments  were  run,  each  1/4  hr,  then  the  composite  accuracy 
would  be: 

Oy  =  1.3/4  =  0.32  arc-sec 

This  is  the  basis  of  the  simplified  one-axis  concept,  measuring  over  a  16-hr 
interval.  Assuming  constant  orientation  relative  to  north,  the  composite 
accuracy  of  a  sequence  of  nor thf indings  would  improve  indefinitely.  Thus,  at  48 
hr,  the  lo  azimuth  error  would  be  0.20  arc-sec,  and  so  on. 


SECTION  5 


ONE-AXIS  SYSTEM  STUDY 


ERROR  ANALYSIS 

The  results  of  the  three-axis  navigator  system  studies  suggest  that  it  is 
pointless  to  measure  rate  around  the  vertical  axis.  We  need  a  different 
approach  and  philosophy.  In  the  horizontal  plane,  the  best  quality  earth-rate 
information  is  conveyed  along  the  east-west  axis.  Earth-rate  measured  around 
the  north-south  axis  is  maximum  for  horizontal  axis  measurements,  and  hence 
varies  only  second-order  with  first-order  azimuthal  drift,  i.e.,  conveys  zero 
azimuthal  drift  information.  Thus,  we  are  led  to  the  concept  of  using  a  single 
gyro,  oriented  east-west.  This  is  the  basis  of  the  one-axis  alternate  system. 

Because  of  the  rather  simple  theoretical  structure,  error  effect  analyses  are 
straightforward,  the  hardest  part  is  deciding  on  a  model  of  the  kinds  of  errors 
to  consider.  For  completeness,  we  assume  that  there  are  random  errors  in 
measuring  the  g,  the  g’,  and  the  e  vectors,  not  necessarily  uncorrelated.  There 
are  seven  variable  errors  to  consider,  namely,  in  el,  e2  and  e3 ,  in  gl  and  g2 , 
and  in  g '1  and  g'2.  g3  is  assumed  to  vary  only  second-order  around  its  nominal 
value  of  1,  and  isn't  measured.  g'3,  the  time  derivative  of  g3,  is  of  course 
also  second  order  and  is  ignored. 

It  is  convenient  to  discuss  errors  in  terms  of  their  main  end  effects.  In  this 
case,  the  obviously  crucial  end  effect  is  the  misdef ini tion  of  the  earth-fixed 
reference  frame,  described  by  the  orthogonal  triad  i,  j,  and  k,  along  east, 
north,  and  upward  vertical.  The  error  end  effect  comprises  small  increments, 
di,  dj ,  and  dk,  respectively  perpendicular  to  i,  j,  and  k.  The  i,  j,  and  k  axes 
are  constructed  from  the  computed  wander-corrected  unit  vector  earth  spin  axis, 
s,  and  from  the  computed  unit  vertical  vector,  v.  For  k,  we  use  v.  Knowing  the 
latitude,  L,  we  express  s  as: 


s  *  j*cos(L)  +  k*sin(L) 


For  j,  we  use  the  k-free  component  of  s,  normalized  to  nnity: 
j  =  (*-k*(k.s))/cos(L) 

(where  ''k.s''  denotes  '’scalar  product  of  vectors  k  and  s' ') 

For  i,  we  simply  vector-crossprodnct  j  and  k: 
i  =  (s  x  k)/cos(L) 

The  error  in  k  is  simply  the  error  in  v,  conveniently  expressed  as 


dk  =  dv  =  i*vl  +  j*v2 

The  error  in  s  is  conveniently  expressed  in  terms  of  east-west  and  meridional 
components,  along  i  and  along  k*cos(L)-j *sin(L) ,  respectively: 

ds  =  i*sl  +  (k*cos(L) -j*sin(L) ) *s2 

The  errors  in  i,  j,  and  k,  in  terms  of  dv  and  ds  are 

di  =  ( (ds  x  k)  +  (s  x  dv))/cos(L) 

dj  =  (ds  -  dv*(k.s)  -  k*(dv.s)  -  k*(k.ds) ) /cos(L) 


These  effects  can  be  represented  by  a  rotational  error  vector,  z: 

di  =  z  x  i 
dj  =  z  X  j 
dk  =  z  z  k 

where: 

z  =  -i*v2  +  j*vl  +  k*(vl*sin(L)-sl) /cos(L) 

The  ref erence- frame  errors  depend  on  the  two  components,  vl  and  v2,  of  the 
vertical  error  and  on  the  east-west  component,  si,  of  the  earth  spin  axis 
error.  Since  the  meridional  component  of  the  earth  spin  is  ignored,  by  relying 
solely  on  the  northr-south  component  of  the  tiltmeters,  s2  plays  no  role  in  the 
error  analysis.  The  rationale  for  so  doing  is  the  assumption  that  the  tiltmeter 
accuracy  will  be  considerably  better  than  the  gyro  accuracy,  so  axy 
inconsistency  between  tiltmeter  readings  and  gyro  readings  (implying  a  false 
value  of  latitude)  is  to  be  resolved  in  favor  of  the  tiltmeters. 

The  error  components  vl ,  v2 ,  and  si  are  simply: 

vl  =  err(a2)  =  err(gl) 
v2  =  err(-al) =  err(g2) 
si  =  -err(dl/R) 

=  sin(L)*err(gl)  -cos(L)*err(a3)  -  (err(el)  +  err(g'2))/R 

Substituting,  the  rotational  error  vector  becomes: 

z  =  -i*err(g2)  +  j*err(gl) 

+k*(err(a3)  +  (err(el)  +  err(g *2) )/(R*cos(L) ) ) 
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If  we  use  independently-generated  values  of  earth  spin  axis  deviation,  or 
multi-site  averaging  of  vertical  axis  warp,  we  dispense  with  the  theodolite 
observations  of  landmarks,  and  the  a3  contribution  to  z  disappears,  giving: 

z  =  -i*err(g2)  +  j*err(gl) 

+k*(err(el)  +  err(g *2) ) /(R*cos(L) ) 

This  says  that  the  tilt  component  of  the  rotational  error  (i.e.,  that  part  which 
is  around  a  horizontal  axis),  is  due  to  the  tilt  measurement  error 
(crossproducted  with  the  unit  vertical),  and  that  the  azimuthal  rotation 
component  is  due  to  east-west  axis  rotational  error,  comprised  of  the  sum  of 
gyro  east-west  axis  rate  measurement  error  and  of  east-west  axis  rate-of-tilt 
measurement  error.  The  foregoing  results  are  completely  in  harmony  with 
intuition.  The  components  of  z  around  horizontal  axes  are  governed  simply  and 
solely  by  the  tiltmeter  errors.  The  component  of  z  arou-J  the  vertical  axis  has 
the  expected  sensitivity  to  latitude.  The  errors  in  gl  and  g2  are  respectively 
east-west  and  north-south  tilt  errors,  and  that  in  g '2  is  the  error  in 
determining  the  rate  of  change  of  north-south  tilt.  The  error  in  el  is  the 
error  in  measuring  east-west  rate. 

GYRO  PERFORMANCE  REQUIREMENTS 

Assume  the  performance  is  to  be  dominated  by  gyro  random  drift  error.  This  is 
justified  by  the  high-quality  state  of  the  art  in  tiltmeters  and  by  the  fact 
that  we  can  remove  the  effect  of  gyro  low  frequency  bias  errors  by  indexing  (as 
well  as  tiltmeter  bias  effects).  Then,  the  azimuthal  component  of  the  error 
vector,  z,  wil  be  equal  to: 

z(az'l)  =  err(el) *secant ( la ti tude ) / ear thra te 


The  rms  error  in  el  will,  under  angle  random  walk  (which  is  rate  white  noise), 
be  equal  to  the  random  walk  coefficient  divided  by  the  square  root  of  the  time 
of  observation,  giving: 

z(az'l)  =  RW*SEC(L)/(R*SQRT(TAO)) 

where: 

RW  =  Random  drift  coefficient,  deg/rt-hr 
L  =  Latitude 
R  =  Earthrate,  deg/hr 
TAU  ■  Time  of  observation,  hr 

Putting  R  =  15.0411°/hr,  and  assuming  a  measurement  time  of  16  hr,  and 
converting  z(az’l)  from  radians  to  arc-seconds: 

z(az'l)  =  20626  5*RW*SEC(L) /(IS.  0411*4) 

=  3428*SEC(L)  *R!f  arc-sec  rms 

Assuming  that  we  want  z(az'l)  -  0.2  arc-sec,  the  random  walk  value  needed  to 
achieve  this  is 

RW  «  0.2*COS(L)/3428 

=  0. 00005 8*C0S(L) deg/ rt-hr 

Thus,  at  the  equator,  we  need  a  0.000058  gyro.  At  latitude  45°,  0.00004  would 
be  needed. 


EXPECTED  SYSTEM  PERFORMANCE 

The  finite  effects  of  bias  instability  and  quantization  are  considered 
explicitly  as  follows.  With  indexed  readings  taken  over  total  time  T,  with  each 
indexed  reading  lasting  time  t,  the  total  xms  error  angular  error  is  given  by 
the  expression: 

SQRT(QS**2/(t*T)  +  RW**2/T  +  3* (BI*t/T) **2) /(R*COS(LT) ) 


where : 


QS  =  ms  quantization/ spillover 
RW  -  Random  walk  coefficient 
BI  =  Low  frequency  bias  instability 
R  =  Earth  rate 
LT  =  Latitude 

T  =  Total  measurement  time 
t  =  Indexation  period 


=  0.024  arc-sec 
*  0. 00004°/ rt-hr 

-  0.0003°/hr 
=  15.041°/hr 

-  45° 

-  16  hr 


The  gyro  parameter  values  are  from  the  experimental  undithered  gyro  data 
reported  in  Section  6.  The  only  assumption  is  that  dithered  operation  can  be 
enhanced  to  give  the  same  random  walk  as  undithered  operation,  using  closed- 
loop  lockin  correction  (CLIC)  to  cancel  out  random  errors  accrued  at  each  dither 
reversal,  as  discussed  in  a  later  subsection  (Gyro  Modifications  Needed). 


Using  the  above  numerical  values,  the  rms  error  in  arc-sec,  as  a  function  of  t, 
is  given  by: 

ms  error  «  SQRT(0.001/t  +  0.038  +  0.4*t«*2) 


This  has  a  minimum  value  of  0.23  arc-sec  at  t  >  0.11  hr.  Therefore,  the  optimum 
indexation  rate  is  about  9  times  per  hour. 


HARDWARE  REQUIREMENTS 


Since  only  the  el  component  of  rate  error  matters,  it  makes  sense  to  minimize  el 
by  using  all  the  time  available  to  take  east-west  rate  readings.  This  suggests 
that  we  need  only  an  east-west  gyro.  An  added  benefit  of  using  the  gyro  along  a 
nominally  null-rate  axis  is  that  gyro  scale-factor  variability  is  not  an  issue. 
Similarly,  tiltmeter  scale  factor  variability  is  removed  from  the  picture  by 
installing  the  tiltmeters  nominally  horizontal. 

Factors  contributing  to  the  error  in  determining  the  el  component  of  the  earth 
spin  vector  are  gyro  bias  error  and  gyro  input  axis  error,  as  follows: 

err(el)  *  dB  +  R*(dV*sin(L)  +  dH*cos(L)) 

where : 

dB  =  Deviation  of  gyro  bias  rate  from  nominal 
dV  =  Vertical  deviation  of  gyro  axis  from  nominal 
dH  =  Horizontal  deviation  of  gyro  axis  from  nominal 

If  the  deviations,  though  unknown,  are  constant  over  the  observation  period,  we 
can  cancel  their  effects  by  taking  indexed  readings,  pointing  first  east,  then 
west,  and  evaluating  the  east  axis  rate  as  the  averaged  east  reading,  minus  the 
averaged  west  reading,  divided  by  2.  Indexing  around  the  vertical  axis  will 
cancel  the  error  contributions  from  dB  and  d V.  Indexing  around  the  north-south 
axis  will  cancel  the  error  contributions  from  dB  and  dH.  If  the  deviations  are 
time-variant,  the  benefits  of  indexing  will  depend  on  the  dynamics  of  the  time 
variation. 

High  frequency  bias  variation  is  exemplified  by  ordinary  gyro  angle  random  walk 
i.e.,  rate  white  noise,  against  which  indexing  has  no  effect.  Hence,  gyro 
output  angle  random  walk  is  a  major  driver  in  this  application.  Laser  gyro 
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input  axis  variation  is  believed  to  be  of  order  1  arc-sec  and  to  vary  very 
slowly,  with  a  time  constant  of  hours,  at  least.  Thus,  indexing  should  greatly 
dilute  the  dV  contribution,  probably  to  a  negligible  level. 

Indexing  introduces  it's  own  set  of  added  errors,  corresponding  to  the  accuracy 
and  repeatability  with  which  the  indexed  orientations  are  achieved.  The 
UL TRADES  commercial  indexing  device  can  index  back  and  forth  through  180°  with 
error  of  <  1/8  arc-sec  and  repeatability  better  than  1/20  arc-sec,  around  the 
index  axis.  The  tilt  error  is  ostensibly  of  order  1  to  2  arc-sec,  and  is  some 
variable  blend  of  random  from  one  index  to  the  next  and/or  quasi-constant. 

Any  random  component  would  be  attenuated  over  N  successive  indexes,  by  a  factor 
of  SQRT(N).  Also,  it  ought  to  be  possible  to  directly  measure  the  tilt  of  the 
ULTRADEX'd  upper  surface  (bearing  the  gyro  and  tiltmeters)  relative  to  the  lower 
main  stationary  surface  on  which  the  ULTRA DEX  rests.  In  principle,  sensitive 
height  gages  could  measure  the  elevation  of  the  upper  surface  over  the  lower,  at 
four  equi-spaced  points  around  the  rim  of  the  upper  surface.  A  trivial 
calculation  would  yield  the  two-axis  relative  tilt. 

GYRO  MODIFICATIONS  NEEDED 

With  only  one  gyro,  it  is  impractical  to  use  kinematics  such  as  are  used  in  the 
three-axis  system.  Three-axis  rate  information  is  a  must  for  talleying  compound 
rotation.  A  one-axis  system  must  keep  the  gyro  in  fixed,  or  highly  constrained, 
orientations  to  avoid  uncomputability  situations.  But,  the  whole  point  of  the 
three-axis  kinematics  is  to  keep  the  gyros  out  of  lockin.  At  very  low  input 
rates,  the  laser  gyro  scale  factor  degrades  to  zero,  resulting  in  output  error. 
So,  some  other  kinematics  must  be  used,  which  both  keeps  the  gyro  out  of  lockin 
and  preserves  computability.  The  standard  way  of  doing  this  is  to  oscillate  the 
gyro  around  its  input  axis,  referred  to  as  dithering  the  gyro.  This  avoids  the 
computability  problems  of  compound  motion,  exemplified  by  coning  errors. 


Dithering  is  conveniently  done  by  mounting  the  gyro  on  a  torsional  spring 
(oriented  around  the  input  axis)  and  driving  the  gyro  in  an  angular  sinusoidal 
motion,  typically  a  few  hundred  arc-sec  peak-to-peak,  at  a  frequency  of  a  few 
hundred  hertz. 

Now,  there  is  another  problem.  The  dither  oscillation  goes  to  zero  rate  twice 
every  dither  cycle,  and  during  these  brief  instants,  lockin  effects  cause  the 
gyro  output  to  accrue  small  increments  of  error.  Their  rms  value  is  a  function 
of  intrinsic  gyro  lockin  rate  and  of  the  time  spanning  each  dither  reversal 
during  which  the  angular  dither  rate  is  less  than  the  lockin  rate.  Under 
ordinary  randomized  dither,  the  errors  accumulate  to  a  random  walk  angular 
output  error.  Typical  dithered  gyro  random  walks  are  0.003.  0.001,  and 
0 .00025°/ rt- hr,  for  the  GG1342,  GG1300,  and  GG1389  Honeywell  RLGs, 
respectively.  As  far  as  we  know,  the  GG1389  has  the  lowest  ordinary-dither 

random  drift  of  any  laser  gyro  in  the  world,  but  its  value  is  still  far  too  high 

to  enable  subarc-sec  system  performance. 

Analytic  modeling  of  the  laser  gyro  has  disclosed  that  the  error  increments 
incurred  across  each  dither  reversal  are,  other  things  being  equal,  sinusoidal 
functions  of  the  nominal  value  of  the  gyro  phase  angle  at  the  dither  peak. 

Under  randomized  dither  (in  which  the  dither  amplitude  varies  slightly  for  each 

successive  dither  peak),  a  random  sequence  of  error  increments  is  generated, 
which  adds  to  a  random  walk.  This  model  has  been  confirmed  by  many  experiments. 

Observed  random  drifts  agree  very  well  with  this  model  (which  is  is  based  on 
approximating  the  sinusoidal  dither  near  the  peak  as  a  parabola,  and  then 
integrating  the  lockin  equation  to  yield  a  pair  of  Fresnel  integrals  to 
represent  the  error  increment).  So,  it  was  a  natural  step  to  pursue  the  model 
further  and  find  ways  to  reduce  the  random  drift.  The  CLIC  invention  was  the 
result.  This  functions  by  controlling  the  dither  in  detail,  in  such  a  way  that 
the  trig  functions  (sine  and  cosine)  of  the  gyro  phase  angle  values  at  the 
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dither  minimizes  sum  to  zero,  rather  than  to  an  unbounded  random  walk.  Then, 
with  the  trig  functions  held  to  zero  sum,  the  drift  error  will  likewise  become 
bounded  and  low. 

Upon  implementing  CLIC,  gyro  random  drifts  did  decrease  by  a  factor  of  up  to 
tenfold.  With  some  ezceptions,  it  appears  that  CLIC  makes  the  gyro  perform  up 
to  the  so-called  quantum  limit,  as  governed  by  the  canonical  uncertainty 
equations  of  quantum  mechanics.  Quantum-limited  performance  is  also  approached 
when  the  gyro  is  continuously  rotated  in  one  direction.  This  has  been  shown 
over  five  years  of  research  and  tests.  The  hoped-for  value  of  0 .00004°/ rt-hr  is 
about  two  times  the  calculated  quantum  limit.  It  seems  reasonable  to  expect 
this  value  from  a  CLIC-enhanced  GG1389,  especially  since  it  has  been  achieved 
under  unidirectional  rotation  testing. 

CLIC  has  a  side  benefit  that  is  also  vital  to  getting  subarc-sec  performance. 

It  provides  enhanced  angular  resolution,  to  within  <  0.01  arc-sec,  instead  of 
the  usual  laser  gyro  pulse  weight  quantization  (1  arc-sec,  for  the  GG1389). 

This  is  done  by  the  CLIC  software  for  the  original  purpose  of  calculating  gyro 
drift  error  discriminant  values.  Without  this  feature,  each  of  the  indexed 
readings  taken  would  be  corrupted  by  1/SQRT($)  gyro  pulses,  i.e.,  0.4  arc-sec 
rms.  Even  with  only  four  readings  per  hour,  this  would  generate  added  error  of 
0.00020°/rt-hr.  This  would  result  in  performance  of  5  times  spec,  instead  of 
being  almost  in  spec.  But,  with  0.01  arc-sec  resolution,  there  would  be 
negligible  added  effect. 

CLIC  has  been  well  developed  and  refined  over  the  last  five  years  and  needs  no 
further  development,  beyond  putting  it  into  a  different  gyro  than  the 
presently-used  GG1342.  This  entails  some  hardware  additions  to  the  gyro  dither 
drive  and  readout,  and  simple  changes  in  the  software — mainly  different  values 
of  constants  for  gyro  scale  factor,  dither  frequency,  drive  scale  factor,  etc. 
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SECTION  6 


EXPERIMENTAL  RESULTS 

UNDITHERED  GYRO  RANDOM  DRIFT  CHARACTERISTICS 

Having  determined  that  a  snbarc-sac  system  would  require  random  drifts  of  less 
than  0. 000045°/ rt-hr,  we  see  that  only  the  Honeywell  GG1389  RLG  could  meet  this 
goal,  and  then  only  when  continuously  rotated  or,  if  dithered,  only  with  CLIC. 
Since  we  don't  have  CLIC  yet  for  the  GG1389,  we  proceeded  to  concentrate  on 
rate-biased  tests.  These  were  done  by  mounting  the  GG1389  on  our  Goertz  Series 
800  Direct-Drive  Rate  Table,  with  the  input  axis  vertical,  and  rotating  at  30 
and  40°/sec. 

The  table  has  bounded-error  in  its  angular  measurement,  of  a  nature  that  is 
readily  dealt  with  by  processing  the  output  with  a  triangular  filter.  This 
drastically  attenuates  the  bounded  spillover/ quantization  error  of  the  table 
readout,  relative  to  the  integrated  white-noise  angle  random  walk,  and  enables 
measurement  of  very  low  random  walks  that  would  otherwise  be  completely  obscured 
by  table  bounded-error  noise. 

The  original  copious  data  stream  obtained  over  2°  or  3°  sectors  of  revolution 
was  compacted  while  being  generated  into  a  much  smaller  number  of 
triangular-f il tered  one-revolution  data  ensembles,  which  were  stored  in  our  Data 
General  computer.  They  were  subsequently  computer-processed  to  generate 
variance  plots,  and  best-fitted  to  a  four-parameter  curve,  comprising  residual 
bounded  error,  angle  random  walk,  bias  instability,  and  rate  random  walk,  as 
defined  by  the  values  of  the  fitting  coefficients. 


Hie  only  thing  that  has  delayed  the  application  of  CLIC  in  systems  is  that  it 
depends  on  a  computerized  prediction  of  how  the  inertial  rotation  input  to  the 
gyro  increments  over  one  dither  cycle  (about  3  msec,  for  the  GG1342) .  It  needs 
to  know  this  to  within  0.2  arc-sec,  in  order  to  generate  the  right  amount  of 
modulation  to  change  the  next  dither  peak  angle  to  a  value  that  helps  to  keep 
the  aforesaid  sums  of  trig  functions  bounded.  So,  if  there's  noise  in  the 
angular  environment  of  order  0.2  arc-sec  uncorrelated  over  3  msec,  CLIC  will 
break  down.  Because  of  this,  CLIC  is  presently  limited  to  benign  environments, 
such  as  in  the  laboratory,  on  a  rigid  mounting. 

The  angular  environment  in  the  situation  being  discussed  here  is  (or  can  easily 
be  made  so),  an  extremely  benign  one.  This  is  a  perfect  application  for  CLIC. 

A  CLIC-enhanced  dithered  laser  gyro  may  be  expected  to  give  random  drift 
performance  equal  to  that  of  elaborately  carouseled  undithered  laser  gyros. 

Now,  with  no  carousel ing,  one  doesn't  need  three  gyros  to  sort  out  the 
kinematics  of  compound  rotations  and  to  compute  platform  attitude.  Rather,  one 
gyro  is  rigidly  attached  in  a  well-defined  relationship  to  the  surface  being 
measured.  A  one-gyro  system  may  not  be  only  equal  to  the  three-axis  system,  but 
even  better. 
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The  results  of  this  testing  procedure  are  shown  in  Figures  10,  11  and  12.  These 
comprise  tests  on  GG1389  gyro  no.  003,  done  on  three  different  days,  at  40°/sec 
with  3°  sector  data,  and  at  30°/sec  with  2°  sector  data.  In  all  these  tests, 
the  best  fit  is  obtained  with  an  angle  random  walk  value  of  0 .000040°/rt-hr.  We 
therefore  infer  this  as  the  official  value  of  undithered  angle  random  walk  for 
this  gyro.  This  is  close  to  meeting  the  requirements  for  the  performance  goal 
of  0.2  arc-sec  rms  measurement  error  over  16  hr,  defined  in  Section  2,  and  is 
even  lower  than  the  value  of  0.000045  assumed  in  Section  4  for  the  Dual-Freedom 
System  Analysis  and  Simulation,  using  GG1389s. 

INP0T  AXIS  STABILITY  TESTS 

Gyro  input  axis  stability  is  as  important  a  factor  as  random  drift  in  a 
subarc-sec  system,  though  of  lesser  relative  importance  in  other  kinds  of 
application.  For  this  reason,  there  had  been  little  formal  investigation  of 
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Figure  10.  GG1389  Random  Drift  Analysis,  17  January  1984 


Figure  12.  GG1389  Rendon  Drift  Analysis,  20  January  1984 


laser  gyro  inpat  axis  stability.  Bat  recently,  this  situation  has  been 
remedied.  We  did  a  series  of  tests  between  March  and  May  1984  on  three  GG1389 
laser  gyros.  These  tests  indicated  a  most  gratifying  degree  of  input  axis 
stability,  relative  to  the  gyro  block  itself,  with  the  worst  component  of  the 
worst  gyro  being  only  +0.04  arc-sec.  This  is  the  intrinsic  limiting  stability. 
The  input  axis  stability  relative  to  the  fixturing  on  which  the  gyros  were 
mounted  ranged  from  0.25  to  0.90  arc-sec,  from  remount  to  remount.  This 
emphasizes  the  need  for  carefully-made  indexing  fixtures  for  the  one-axis 
subarc-sec  demonstrator  unit. 

The  experimental  configuration  was  as  follows:  the  gyro  was  mounted  on  a 
fixture  with  its  input  axis  nominally  in  the  horizontal  plane.  The  fixture 
itself  was  mounted  on  the  Goertz  table,  with  the  latter  having  rotation  around  a 
nominally  vertical  axis.  The  gyro  was  mounted  on  the  fixture  in  four  positions, 
indexed  in  90°  increments  of  rotation  around  the  nominal  gyro  axis.  The  four 
positions  are  denoted  NA  up,  NA  down,  LA  up,  and  LA  down.  LA  and  NA  denote  two 
gyro  body  axes  perpendicular  to  the  nominal  gyro  axis  (see  Figures  13  and  14). 
The  gyro  was  operateu  with  dither  on,  and  with  angular  input  generated  by 
rotating  the  Goertz  table  through  many  revolutions,  both  clockwise  and 
counterclockwise. 

Nominally,  with  the  gyro  axis  perfectly  orthogonal  to  the  Goertz  table  rotation 
axis,  the  only  gyro  output  will  be  gyro  bias  and  random  walk,  plus  a  small 
contribution  due  to  earth  rate,  seen  as  a  sinusoidal  signal,  at  the  frequency  of 
rotation  of  the  Goertz  table,  with  an  amplitude  equal  to  horizontal-plane  earth 
rate  (earth  rate  times  cosine  of  latitude).  But,  with  data  taken  over  an 
integral  number  of  Goertz  table  rotations,  the  latter  sums  to  zero.  (This 
assumes  a  constant  rate  of  rotation  of  the  Goertz  table,  which  is  almost  exactly 
the  case.) 


So,  the  output  will  comprise  bias,  pins  random  drift,  pins  a  contribution  due  to 
the  gyro  axis  unit  vector  having  a  small  component  along  the  Goertz  table 
rotation  axis  unit  vector.  Using  small  angle  approximations,  the  resulting 
output  per  Goertz  table  revolution  will  equal  the  gyro  scale  factor,  in 
counts/revolution,  times  the  component  of  the  gyro  axis  along  the  Goertz  axis, 
times  an  adjustment  for  the  earth  rate  component  along  the  Goertz  axis. 

From  the  results  of  the  eight  runs,  we  can  determine  gyro  bias,  two  components 
of  non-or thogonali ty  of  the  gyro  sensing  axis  to  the  LA-NA  plane,  and 
non- orthogonal ity  between  the  normal  to  the  LA-NA  plane  and  the  Goertz  table 
rotation  axis.  With  this  factor-of-two  overdeterminacy,  we  can  also  form 
estimates  of  the  stability  of  the  above  four  variables.  Hie  computations  are 
based  on  the  following  model  of  gyro  input  axis  deviation. 

Gyro  axis  unit  vector,  V,  =  I  +  perturbation,  D.  The  values  of  D  in  the  4  gyro 
orientations  are 

D(NAU)  =  J*GL+K* ( MK+GN) 

D(NAD)  =  -J  *GL+K* (MK-GN) 

D(LAU)  «  -J*GNH-K*(MK+GL) 

D(LAD)  =  J  *GN*-K*  (MK-GL) 

where : 

I  *  nominal  V,  perpendicular  to  K 
J  =  in-plane  axis,  perpendicular  to  I  and  K 
K  =  Goertz  rotation  axis 
GL  =  Component  of  D  along  LA  body  axis 
GN  *  Component  of  D  along  NA  body  axis 
MK  *  Component  of  tilt  of  indexing  axis  along  K 


The  one- revolution  output  that  would  be  seen  along  the  Goertz  axis  at  rate  R  is 
CR  =  SF*(EV+R)/|R| 

where: 

SF  =  Gyro  scale  factor,  in  counts/revolution 
EV  3  Component  of  earth  rate  along  the  Goertz  axis 
R  =  Rate  of  rotation  around  Goertz  axis 

The  gyro  output  over  one  revolution  will  equal  CR  times  the  component  of  V  along 
the  Goertz  axis,  plus  gyro  bias  rate  times  the  period  of  revolution  times  scale 
factor : 

GY  *  SF*BG*TAD/2*PI  +  DK*CR 

where: 

GY  =  One-rev  gyro  output,  in  counts 
BG  =  Bias  rate  of  gyro,  in  radians  per  second 
TAU  =  Period  of  revolution  around  Goertz  axis,  in  seconds 
DK  3  Component  of  D  along  Goertz  axis 
PI  3  3.1415937022 . 

The  foregoing  is  applied  as  follows:  the  method  will  be  exemplified  by  working 
out  the  set  of  stability  values  for  GG1389  gyro  no.  Oil,  using  the  measurements 
taken  for  that  gyro. 

Equating  the  expressions  for  the  eight  outputs  to  the  values  actually  measured, 
we  get  the  following  eight  equations. 


GY(NAD.CCW) 

s: 

BT  + 

SF*( 1+EPS) *(MK+GN) 

= 

110.159 

counts 

GY(NAU,  CW) 

s: 

BT  - 

SF*(1-EPS)*(MK+GN) 

= 

-110.600 

counts 

GY ( NAD,  CCW) 

= 

BT  + 

SF*( 1+EPS) •(MK-GN) 

- 

-19.560 

counts 

GY(NAD, CW) 

= 

BT  - 

SF* ( 1-EPS )* (MK-GN) 

20.080 

counts 

GY(LAD.CCW) 

= 

BT  + 

SF*< 1+EPS) *<MK+GL) 

- 

33.980 

counts 

GY(LAU,  CW) 

= 

BT  - 

SF*< 1-EPS) *(MI+GL) 

= 

-33.580 

counts 

GY(LAD.  CCW) 

= 

BT  + 

SF*( 1+EPS) *(MK-GL) 

- 

82.220 

counts 

GY (LAD, CW) 

= 

BT  - 

SF* ( 1-EPS) *(MK-GL) 

> 

-81.700 

counts 

where: 

BT  *  SF*BG*TAD/ 2*PI 
EPS  =  EV/ |R| 

Using  Equations  6-1  and  6-3: 

BT+SF* ( 1+EPS) *MK  =  ( 110 .159-19.560) /2  =  45.2995  counts 


Using  Equations  6-2  and  6-4: 


BT-SF*(  1-EPS) *MK  =  -45.2600  counts 


whereby 


MK  -  45 .27 98/ SF 


and 


BT  -  0.0198-45. 27 98*EPS 


BT+SF*(1+EPS)*MK  =  58.1000  count* 

Equations  6-6  and  6-8 

BT-SF*(1-EPS)*MK  «  -57.6400  count* 

whereby : 

MK  =  57. 87 00/ SF 

and 

BT  =  0. 2300-57. 8700*EPS 

Values  of  the  constants  are 

SF  =  1296000  counts/revolution 
R  =  3  0°  /  s  e  c 

EV  =  10.63564°/hr  at  latitude  45°  north 

Hence,  EPS  =  10.63564/(30*3600)  =  0.00009848 

Therefore: 

MK(1234)  =  45.2798/1296000  rad  -  7.20  arc-sec 
MX( 5678)  -  57.8700/1296000  rad  =  8,10  arc-sec 


i.e.,  MX  =  7.65  +0.45  arc-sec. 


BT(1234)  =  0. 0153*180*3600/ 1296000*PI  -  0.0024  arc-sec 
BT( 5678)  »  0 . 2300*1 80*3600/12 9600 0*PI  -  0.0370  arc-sec 

i.e.,  BT  *  0.020  +0.017  arc-sec. 

Furthermore,  subtracting  Equation  6-3  from  Equation  6-1. 

GN( 13)  =  64.859/ <SF*(1+EPS)>  -10.322  arc-sec 
and  subtracting  Equation  6-4  from  Equation  6-2, 

GN(24)  -  65.320/ (SF*(1-EPS))  =  10.397  arc-sec 

i.e.,  GN  -  10.360  +0.037  arc-sec. 

Similarly,  subtracting  Equation  6-7  from  6-5, 

GL( 57)  -  -24 .12/ ( SF*( 1+EPS) )  =  -3.838  arc-sec 

and,  subtracting  Equation  6-8  from  6-6, 

GL( 6 8)  -  -24.06/ ( SF*(1-EPS) )  -  -3.830  arc-sbc 

i.e.,  GL  -  -3.834  +0.004  arc-sec 


The  results  for  no.  Oil  are  suaimarized  in  Table  2. 


TABLE  2.  SHORT-TERM  AXIS  STABILITIES  OF  GG1389  NO.  Oil 


Term 

Value  in 
Arc-Sec 

Variation 
in  Arc-Sec 

ME 

7.65 

±0.45 

BT 

0.020 

±0.017 

GN 

10.360 

±0.037 

GL 

-3.834 

±0.004 

We  saw  a  qualitative  repeat  of  these  results  on  numbers  009  and  004.  That  is, 
very  low  variability  in  three  terms,  and  relatively  high  variability  in  the  tilt 
of  the  fixture  indexing  axis,  ME.  This  further  strengthens  the  conjecture  that 
there  is  either  aniso-sag  of  the  gyro  dither  post  under  the  cantilevered  weight 
load  of  the  gyro  (as  between  the  set  of  NA  up  and  down  measurements  and  the  set 
of  LA  up  and  down  measurements),  or  else  aniso-f lexure  of  the  gyro  base  plate, 
or  even  (less  plausibly)  aperture-connected  bias  shift  due  to  aniso-sag  effects 
on  the  gyro  mirrors,  and  hence  on  the  location  of  the  laser  beam  path  triangle. 

After  the  above  tests,  a  series  of  tests  were  run  over  a  three-day  period, 
comprising  CCW  and  CW  runs  in  the  LA  down  orientation  only.  The  results  are 
presented  in  Table  3. 


TABLE  3.  MOLTIDAY  DATA  OF  GG1389  NO.  Oil 


CCW 

CW 

Sum/ 2 

Diff/4*PI 

90.819 

-89.920 

mam 

92.400 

-91.401 

1 

MM 

92.440 

-91.520 

1 

14.64 

The  daily  readings,  when  differenced,  yield  the  sets  of  day-to-day  differences 
listed  in  Table  4. 


TABLE  4.  MDLTIDAY  STABILITIES  OF  GG1389  NO.  Oil 


Delta- Sum/ 2 


Del ta-Dif  f/ 4*PI 


Day  2  -  Day  1: 
Day  3  -  Day  2: 


-0.04 


Diff/ 4*PI  represents  the  value  of  MX-GL,  while  Sun/2  represents  BT.  The  above 
data  suggests  that  both  MX  (to  do  with  aniso-sag,  etc.),  and  the  gyro  i  nsic 

bias  were  rather  stable  over  a  three-day  period,  more  so  than  was  seen  in  gyros 
1-009  and  1-004. 

The  counts/rev  data  for  GG1389  gyro  no.  009  are  listed  in  Table  5. 

TABLE  5.  MULTI-ORIENTATION  DATA  FOR  GG1389  NO.  009 


NA  up  CCW 
NA  up  CW 
NA  down  CCW 
NA  down  CW 
LA  up  CCW 
LA  up  CW 
LA  down  CCW 
LA  down  CW 


316.329 

-317.700 

298.698 

-300.060 

254.500 

-255.779 

366.423 

-367.962 


From  these  data,  we  computed  the  set  of  stabilities  listed  is  Table  6 


TABLE  6.  SHORT-TERM  AXIS  STABILITIES  OF  GG1389  NO.  009 


Term 

Value  in 
Arc-Sec 

Variation 
in  Arc-Sec 

MK 

49.3 

±0.25 

BT 

-1.018 

+0.007 

GN 

1.400 

+0.005 

GL 

-8.92 

+0.02 

Then,  multiday  stability  data  were  taken,  as  presented  in  Table  7. 


TABLE  7.  MULTIDAY  DATA  OF  GG1389  NO.  009 


ccw 

CW 

Sum/ 2 

Diff/4*PI 

Day  1 

366.423 

-367.962 

-0.77 

58.44 

Day  3 

364.621 

-366.379 

-0.88 

58.17 

378.859 

-380.403 

-0.77 

60.42 

Day  7 

378.561 

-380.282 

-0.86 

60.39 

Plausibly,  the  multiday  variation  of  +.05  arc-sec  in  the  half  sum  suggests  a 
bias  instability  of  order  0.00008°/hr,  The  small  changes  in  the  differences 
between  Days  1  and  3  and  between  days  6  and  7  (0.27  and  0.03  arc-sec, 
respectively)  are  compatible  with  the  one-day  data  stabilities  listed  in  Table 
6.  The  2.2  arc-sec  jump  between  Days  3  and  6  may  be  due  to  some  shift  in 


the  mounting  hardware,  associated  with  shutting  down  for  the  weekend.  The 
GG1389  no.  004  was  also  tested,  giving  the  one-rev  data  listed  in  Table  8. 


TABLE  8.  MULTI-ORIENTATION  DATA  FOR  GG1389  NO.  004 


NA 

up 

CCW 

35.280 

NA 

up 

CW 

-36.000 

NA 

down 

CCW 

198.740 

NA 

down 

CW 

-199.921 

LA 

up 

CCW 

20.120 

LA 

up 

CW 

-21.420 

LA 

down 

CCW 

236.979 

LA 

down 

CW 

-237.980 

The  multiday  data  obtained  are  presented  in  Table  9. 


TABLE  9.  MULTIDAY  DATA  OF  GG1389  NO.  004 


CCW 

CW 

Sum/ 2 

Diff/4*PI 

Day  1 

240.040 

-241.021 

■PVPW 

38.28 

Day  2 

237.501 

-238.478 

BaIw 

37.88 

Day  3 

234.681 

-236.039 

37.46 

From  the  data  in  Tables  8  and  9,  we  infer  the  stabilities  for  GG1389  no.  004 
that  are  listed  in  Tables  10  and  11. 


TABLE  10.  SHORT-TERM  AXIS  STABILITIES  OF  GG1389  NO.  004 


Term 

Value  in 
Arc-Sec 

Variation 
in  Arc-Sec 

MX 

7.65 

±0.45 

BT 

0.020 

+0.017 

GN 

10.360 

±0.037 

GL 

-3.834 

±0.004 

TABLE  11.  MULTIDAY  STABILITIES  OF  GG1389  NO.  004 


Delta  (Sum/ 2) 

Delta  (Dif f/4*PI) 

Day  2  -  Day  1 

0.002 

-0.40 

Day  3  -  Day  2 

-0.190 

-0.42 

TILTMETER  CHARACTERISTICS 

Since,  in  a  one-axis  system,  the  two  tiltmeters  could  do  ''2/3  of  the  job'' 

(tilt  around  north  and  east),  it  is  fitting  to  discuss  tiltmeter 
characteristics.  Figure  15  shows  ( inf erential ly)  the  stability  of  the  simple 
bubble-level  tiltmeter,  as  it  existed  in  1963,  over  a  one-year  period.  TVro 
bubble  levels  were  placed  along  east-west  and  along  north-south,  for  a  total  of 
four.  For  clarity,  the  vertical  scales  were  staggered,  so  that  each  of  the  four 
outputs  could  be  clearly  seen.  The  two  pairs  track  each  other  very  well.  Apart 
from  some  exceptional  periods  (notably,  around  June  for  the  east-west  pair)  they 
seem  to  track  within  a  fraction  of  an  arc-sec,  perhaps  within  0.1  arc-sec. 
Whoever  prepared  this  graph  would  have  been  a  hero  had  he/she  thought  to  plot 
the  differences  between  the  tiltmeter  pairs!  Unfortunately,  the  data  records 
are  long  gone. 
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These  were  1963  state-of-the-art  bobble  tiltmeters.  There  have  been  many 
refinements  since  then,  notably  in  the  area  of  digitized  readout  mechanizations 
with  accuracies  quoted  at  0.05  are-sec  rms.  Thus,  it  is  reasonable  to  assume 
that  tiltmeter  errors  will  RSS  negligibly  with  gyro  errors. 


SECTION  7 


CONCLUSIONS  AND  RECOMMENDATIONS 

CONCLUSIONS 

1.  It  looks  marginally  feasible  to  measure  crustal  rotation  around  the 
vertical  to  an  accuracy  of  0.2  arc-sec  rms,  over  a  period  of  16  hr,  using 
either  a  three-axis  naviga tor- type  system,  comprising  three  GG1389  laser 
gyros,  plus  three  accelerometers  with  special  two-degree  of  freedom 
carousel ingj  or  a  one-axis  system,  comprising  one  dithered  CLIC-enhanced 
GG1389  laser  gyro,  indexed  between  east  and  west,  plus  two  tiltmeters. 

2.  It  is  impracticable  to  distinguish  between  crustal  rotation  displacement 
around  the  vertical  and  EPA  deviation  along  east-west  using  self-contained 
measurements  from  a  single  apparatus  at  one  site.  Either  a  central  station 
detezmi nation  of  earth  spin  axis  deviation  must  be  made  and  applied  to  the 
total  apparatus  reading,  or  possibly  multisite  averaging  may  be  used. 

3.  The  vertical  axis  component  of  the  three-axis  system  angle  sensing  is 
relatively  useless,  because  it  is  random-walk-driven  without  any 
possibility  of  correction  from  accelerometer  ''tilt’’  readings.  Its  error 
is  bounded  only  by  conflicting  rate  information  generated  along  the  east 
axis. 

4.  The  two  tilt  components  of  earth  crustal  rotation  (around  north  and  east) 
can  be  readily  and  adequately  measured  by  two  tiltmeters.  Gyros  are  needed 
only  to  add  information  around  the  vertical. 

5.  The  best  deployment  of  gyro  power  is  along  east-west. 
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6.  The  one-axis  system  would  probably  work  at  least  as  well  as  the  three-axis 
system,  assuming  that  accurate  indexing  can  be  done. 

RECOMMENDATIONS 

Based  on  the  results  of  this  study  and  on  our  perceptions  of  AFGL's  needs,  we 
recommend  that  AFSL  initiate  a  program  to  design,  build,  and  test  a 
transportable  one-axis  demonstrator .  As  well  as  showing  (or  disproving) 
feasibility,  the  demonstrator  prototype  will  indicate  further  areas  to  be 
studied,  refinements  needed  in  the  kinematics,  added  error  sources,  etc.  It 
will  be  the  first  transportable  subarc-second  crustal  deviation  measurement 
device  to  be  built,  so  all  the  data  it  generates  will  be  of  the  utmost  interest 
in  guiding  future  development. 

Recommended  Demonstrator  Design 

The  proposed  demonstrator  is  a  single-axis  indexed  system  with  two  tiltmeters 
and  a  single  GG1389  RLG.  The  sensors  and  their  electronics  are  fixed  to  a 
platform  that  is  mounted  on  a  computer-driven  UL1RADEX  precision  indexing 
table.  The  outputs  of  all  of  the  sensors  will  be  continuously  monitored  via  a 
dedicated  processor  that  will  perform  the  algorithm  calculation  and  output  the 
results  to  a  graphics  display,  thus  allowing  continuous  monitoring  of  the 
performance  of  the  system. 


A  block  diagram  of  the  demonstrator  hardware  is  shown  in  Figure  16.  The  laser 
gyro  is  mounted  with  its  input  axis  nominally  horizontal,  in  the  plane  of  the 
plate  surface,  along  a  direction  that  will  be  nominally  east-west  when 
installed.  Tmo  tiltmeters  are  mounted  to  sense  tilts  around  nominal  east-west, 
and  around  nominal  north-south.  A  block  schematic  of  the  electronics  is  given 
in  Figure  17. 


PATH  LENGTH  CONTROL 


Figure  17.  Block  Schematic  of  Demonstrator  Electronics 
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The  ULTRADEX  indexing  table  is  the  platform  to  which  all  the  sensors  and  their 
electronics  are  fixed.  This  table  is  computer-precision  controlled,  with 
position  error  of  <  1/8  arc-sec,  and  repeatability  better  than  1/20  arc-sec, 
around  the  index  axis.  The  rotation  axis  instability  is  nominally  1  to  2 
arc-sec  random.  This  latter  component  is  attenuated  over  N  successive  indexes, 
by  a  factor  of  SQRT(N).  If  this  instability  becomes  a  concern,  it  would  be 
possible  to  measure  the  relative  tilt  angle  between  the  top  and  bottom  surfaces 
of  the  platform  and  explicitly  eliminate  the  effect  of  rotation  axis 
instability.  This  would  entail  using  added  high-stability  tiltmeters.  The 
present  system  design  does  not  include  extra  tiltmeters,  but  has  expansion 
capability  to  add  them  at  a  later  date.  In  addition,  the  system  will  allow 
changing  to  different  RLGs  with  minimal  modification,  thereby  facilitating 
4  system  performance  upgrades  as  technology  produces  even  better  RLGs. 

The  computation  algorithm  will  be  fed  by  the  outputs  of  the  sensors,  and  its 
cycle  will  extend  over  each  period  between  successive  indexations.  The  results 
of  successive  cycles  will  be  appropriately  filtered  to  yield  a  progressively 
improved  estimate  of  1)  crustal  tilt  around  some  horizontal  plane  axis,  and  2) 
combined  east-west  deviation  of  the  polar  axis  and  crustal  rotation  around  the 
vertical.  The  control  and  algorithm  processor  planned  for  use  is  the  Hewlett- 
Packard  HP9826.  This  machine  has  been  chosen  for  its  interface  capability  to 
physical  sensors  and  for  its  straightforward  high-order  language  operating 
system.  These  factors  will  allow  the  control  and  algorithm  software  to  be 
written  in  a  concise  and  efficient  manner. 

9 

Demonstrator  Features 


Automated  Operation — Fully  automatic  data  acquisition  and  control,  allowing 
unattended  operation  for  extended  periods  of  time. 


On-Line  Calculation — Real-time  calculation  of  combined  polar  axis  deviation  and 
earth  crustal  rotation  around  the  vertical,  from  the  initial-condition  state. 

Portability — Rugged  design,  combined  with  reasonable  size  and  weight,  to  allow 
easy  transport  to  other  laboratory  sites  for  test  and  demonstration. 

Reliability — RLG  life  and  reliability  have  been  clearly  demonstrated  in 
commercial  aircraft  revenue  service,  with  M7BF  in  excess  of  S0000  hr. 
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APPENDIX  A 


ANALYSIS  OF  DUAL- FREEDOM  ROTATIONAL  STRATEGY 


ERROR  MODEL 

The  following  error  model  is  assumed. 

ft 

*  -  -C  5»  (A-l) 

g 

6V  =  C  5a  -  ^xa  (A-2) 

—  a 

where : 

¥  =  attitude  error  vector  of  orthogonal  ISA  system  axes  relative  to  an 
earth-fixed  frame 
6V  =  velocity  error  vector 
a.  =  nongravitation  acceleration  vector  of  ISA 
Cg  =  transformation  matrix  from  gyro  input  axes  to  local-vertical 
reference  frame 

Ca  =  transformation  matrix  from  orthogonal  accelerometer  input  axes  to 
local-vertical  reference  frame 
5o>  =  gyro  output  error  vector 
5j»  =  accelerometer  output  error  vector 

The  question  is,  how  must  t  sensor  assembly  be  rotated  to  cause  attitude  and 
velocity  error  to  have  a  mea  alue  of  zero. 

The  matrix  Cg  can  be  expressed  as  the  product  of  the  matrix  A  (which  relates 
the  gyro  input  axes  to  the  system  frame),  and  the  matrix  C  (which  relates  the 
ISA  frame  to  the  local-vertical  reference  axes). 

Cg  -  CA  (A-3) 
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The  total  error  in  the  gyro  output  (6(g)  can  be  decomposed  in  terms  of  bias, 
scale  factor,  and  misalignment  errors  as 


6<i>  =  6  A  ^u>  .  ,  +  K  at  ,,  +  b 
-  -r/i  g“g/i  ~g 

where: 

Ujc/i  =  inertial  angular  velocity  of  the  ISA 
Wg /£  *  inertial  angular  velocity  measured  by  the  gyro  triad 
bg  =  gyro  bias  error 

5A~1  =  error  in  transformation  matrix  from  ISA  axes  to  gyro  axes 
Kg  =  gyro  scale  factor  error  matrix 

Putting  Equations  A-3  and  A-4  into  Equation  A-l  yields: 

V  =  CA(6A_1(i»  +  K  <o  .  +  b  ) 

-r/i  g-g/i  ~g 

=  CA(8A_1o»  . .  +  K  A_1w  ..  +  b  ) 

~r/ 1  g  -r / i  ~g 

=  CA[(6A_1  +  K  A_1)  u  ,  .  +  b  ] 
g  — r  /  i  — g 


The  effect  of  gyro  input  axis  misalignment  on  A  *  can  be  expressed  as 


1 


i=l 


8. 

x 


where  the  8^'s  are  the  gyro  misalignment  angles.  The  gyro  scale  factor 
matrix.  Kg,  can  also  be  explicitly  defined  as 

K  =  diag  (k  ,  k  „  k  ,) 
g  gl  g2  g3 


(A-4) 


(A-5) 


(A-6) 


error 


(A-7) 


where  the  kg^s  are  the  gyro  scale  factor  errors  and  '  *diag"  denotes  a 
diagonal  matrix.  Substituting  Equations  A-6  and  A-7  into  Equation  A-S  gives: 


V  =  CA 


an?- 


*  4U‘  (lgl  lg2  V* 


-1)  u  . .  +  b 

J  -r/i  “8 


(A- 8) 


The  accelerometer  output  error  is  expressed  as 

6a  =  K  a  +  b 

—  a  “  ~a 


( A-9) 


where : 

a  =  nongr avitational  acceleration  expressed  with  components  in  the  system 
frame 

ba=  accelerometer  bias 


and  Ka  is  a  matrix  that  accounts  for  both  misalignment  and  scale  factor  error 
(the  diagonal  entries  being  scale  factor  errors  and  the  off-diagonal  entries 
being  misalignment  errors).  Since  the  accelerometer  input  axes  nominally 
coincide  with  the  ISA  system  axes,  the  matrix  Ca  in  Equation  A-2  can  be 
replaced  by  C,  which  allows  Equation  A-2  to  be  written  as 

6V  =  C[K  a  +  b  ]  -  f  x  a  (A-10) 

a  a 


CONDITIONS  FOR  COMMUTATION  OF  GYRO  ERRORS 


Referring  to  Figure  A-l,  we  see  that  the  angular  rotation  vector  of  the  ISA  can 
be  expressed  as 


“r/i  = 


an  - 
~x 


(A-ll) 
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a.  Inertial  frame  b.  ISA  reference  frame  c.  Outer  gimbal  reference  frame 
Figure  Arl .  Definition  of  Coordinate  Frames 

where : 

a  =  outer  gimbal  inertial  rotation  rate 
P  =  inner  gimbal  inertial  rotation  rate 

and  the  inertial  rotational  rate  about  the  axis  perpendicular  to  the  two  gimbal 
axes  is  taken  to  be  zero.  The  unit  vector  uz  belongs  to  the  triad  (uz,  Uy.  uz) 

attached  to  the  outer  gimbal.  and  0*3  belongs  to  the  triad  (urj ,  ur2#  ur3>  that 
defines  the  ISA  axes.  If  these  two  triads  are  assumed  coincident  at  t^O,  the 
subsequent  transformation  between  them  is  defined  by: 


which  allows  the  ISA  rotation  rate  to  be  expressed  as 


id  -  oCpu, 
-r/ 1  1 


aSpu^ 


(/ 


If  the  inertial  reference  triad  (u^j,  u*2»  0^3)  is  also  taken  as  being 
coincident  with  the  other  two  triads  at  t=0 ,  the  following  is  true. 


Which,  with  Equation  A-12,  leads  to  the  defining  relationship  for  C: 


i 

ai 

cp 

sp 

0 

r 

*1 

r 

5i 

i 

*2 

= 

-CaSp 

CaCp 

-Sa 

r 

% 

=  C 
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*3 

-SaSp 

SaCp 
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Using  Equations  A-13  and  A-14  the  attitude  error  propagation  equation  (A-8) 
can  be  expressed  as 


where  &g  is  the  gyro  error  vector: 


e 

“g 
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and  the  elements  of  the  matrix  H  are  given  by: 


h.  . 

ij 


m 


k=l 


(A-16) 


(A-17) 


in  which  f^(t)  is  the  l  th  of  M  time-dependent  functions,  and  the  aj.'  s  are 
constants  that  are  functions  only  of  the  elements  of  the  transformation 
matrix  A. 


If  all  sensor  errors  are  assumed  constant  over  time  T,  then  Equation  A-15  can 
be  directly  integrated: 


V(T) 


T 

/  H(t)dt 
o 


e. 


(A-18) 


Furthermore,  for  ^(T)  to  equal  zero: 

o/Tfk(t)dt  *  0  k  =  1 ,  2,  - M 


( A-19) 


The  integrals  of  the  f^(t)  that  appear  in  the  solution  of  the  attitude  error 
vector  ¥  are  defined  in  Table  A-l,  which  shows  that  all  of  the  functions  to  be 
integrated  depend  only  on  the  gimbal  angles  and  rates.  If  the  ISA  inertial 
rotation  rates  about  the  two  gimbal  axes  are  constant,  all  of  the  integrals  can 
be  evaluated  explicitly.  As  an  example,  one  of  the  Table  A-l  integrals  is 
evaluated  as  follows. 


Originating  from 

Originating  fro«  Cyro  Scale  factor  and  6yr0  BUs  trrors 

lllsallgiwent  trrors 


TABLE  A-l .  INTEGRAL  FUNCTIONS  FOR  ATTITUDE  ERROR  PROPAGATION 


Jf  k(t)dt 

Acceptable  Ratio  of 
fiinhal  Rates 

J Co*8dt 

All  n 

f SlnBdt 

All  n 

^Slnodt 

All  n 

Jcosadt 

All  n 

JsinaSInBdt 

n  j»  1 

J CosaSinSdt 

All  n 

JsinoCosSdt 

All  n 

j CosaCosadt 

n  t  1 

J  4S1n2Bdt 

All  n 

JiSInaSInBdt 

n  t  \ 

Ji CosoS InBdt 

All  n 

JaSInaCosadt 

All  n 

J’aCOSaCosadt 

All  n 

Ji  S1nZadt 

No  n 

JiCos2adt 

No  n 

j'iS1naSin2adt 

All  n 

J’iCosaS1n2adt 

All  n 

JaS1naCos28dt 

All  n 

J  aCoiaCos2sdt 

All  n 

JiS1n,$1n28dt 

All  n 

JiC0SaS1n28dt 

All  n 

•  The  ratio  (n)  of  Inner  to  outer 
glrabal  rates  that  produces  a 
zero  value  of  the  Integral  over 
the  period  of  the  outer  gimbal 
motion 
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TABLE  A-l.  INTEGRAL  FUNCTIONS  FOR  ATTITUDE  ERROR  PROPAGATION  (concluded) 


/V* 

Acceptable  Ratio  of 

Gimbal  Rates 

JaCos8dt 

All  n 

JsSInsdt 

All  n 

v» 

U 

o 

01  u 
—  u 

J'$S  Inadt 

All  n 

«  Ui 

o 

4/1  ■*-» 

c 

o  0> 

b  s 

f SCosadt 

All  n 

o’  Ol 

8  i 

in 

JsSInaSInedt 

n  /  1 

as 

C  T3 

^  c 
«-* 

J"  sCosaSInsdt 

All  n 

c  u 
•—  o 

O  4-» 

u 

i;  * 

J*8  SlnaCossdt 

All  n 

o  u. 

J"  sCosaCossdt 

n  f  1 

T*  2 

/  aCosaSin  fldt 


=  a  /^osatSin^ptdt 


=  2  [o/Tc°S“ 

=  j  |j5inatj 


atdt  -  /  CosatCos2ptdt 

o 


tptdtj 


a 

4 


Sin(a  -2B)  t  Sin(a  +  23)  t 
(i  -  2f3)  (a  +  2p) 


Now,  assume  that  the  inner  gimbal  rate  is  n  times  the  outer  gimbal  rate,  and 
that  the  integration  interval  (T)  is  equal  to  the  outer  gimbal  period  of 
revolution.  Then,  the  integral  can  be  evaluated  as 
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f  -  j  Sid.  -  }  [**%$*  *  “f4] 

=  0  for  all  integer  N 

All  except  two  integrals  in  Table  A-l  reduce  to  zero  over  one  outer  gimbal 
period,  for  integer  N.  The  two  exceptions  are 

aS in  pdt  and  /A  aces'4 pdt 

o  o 

both  of  which  give  aT/2  over  the  period  T.  But,  by  reversing  the  outer  gimb&l 
rate  after  each  successive  revolution  of  360°,  these  two  give  zero  net  over 
intervals  of  2T.  The  periodic  reversal  of  a  has  another  benefit  as  well;  it 
causes  every  integral  that  has  a  as  a  factor  to  have  a  mean  value  of  zero 
(providing  that  any  restriction  on  the  ratio  of  gimbal  rates  is  observed).  The 
result  is  that  attitude  error  will  tend  to  be  more  nearly  zero-mean,  which  in 
turn  produces  a  lower  growth  of  velocity  error  through  the  V  X  a,  term. 

CONDITIONS  FOR  COMMUTATION  OF  ACCELEROMETER  ERRORS 


The  velocity  error  expression  is 

6V  =  We 
_  — * 


( A-20) 


where  is  a  vector  formed  from  the  accelerometer  errors  (bias,  scale  factor, 
misalignment),  and  the  elements  of  the  matrix  W  take  the  general  form: 


w .  . 

tj 


bkgk(t) 


(A-21) 
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in  which  gfc(t)  is  the  t  th  of  N  tine-dependent  functions,  end  the  b^'s  ere  e 
function  of  the  conponents  of  nongrevitetionel  ecceleretion,  a. 


If  ell  accelerometer  errors  ere  teken  to  be  constent  over  en  integretion 
interval  T  (the  period  of  the  outer  giabel  notion),  end  if  the  conponents  of 
nongrevitetionel  ecceleretion  ere  elso  essentially  constent  over  T,  Equation 
A-20  hes  the  solution: 


5V(T)  = 


/TW(t)dt 


e 

a 


(A-22) 


which  produces  a  zero  value  if: 

/Tg. (t)dt  =0  k  -  1,  2,  ...N 

O  K 


The  integral  functions  that  appear  in  the  solution  of  the  velocity  error 
equations  ere  defined  in  Table  A-2,  some  of  which  ere  the  sene  as  those 
appearing  in  the  attitude  error  equations.  As  can  be  seen  from  Table  A-2,  ell 
but  eight  of  the  integrals  can  be  made  to  produce  a  zero  value  at  the  end  of 
each  complete  revolution  of  the  outer  gimbal.  However,  since  none  of  the 
integrals  associated  with  wcelerometer  errors  involve  the  outer  gimbal  rate, 
reversal  of  the  outer  gimbal  motion  will  not  produce  a  zero  value  at  the  end  of 
2T.  Therefore,  the  contribution  to  velocity  error  represented  by  these  eight 
terms  will  be  generally  monotonic.  The  accelerometer  errors  associated  with 
these  integrals  are  scale  factor  and  misalignment  errors,  and  the  integrals 
are  associated  with  error  in  the  vertical  direction  (i.e.,  in  the  direction  of 
the  gravity  vector).  Furthermore,  the  integrals: 

/Sin^aCospdt  /C os^aCospdt  /Sin2aSinPdt 


will  not  produce  a  zero  value  over  T  if  the  ratio  of  inner  gimbal  to  outer 
gimbal  rate  (n)  is  2.  Therefore,  an  n  value  of  3  or  higher  should  be  used. 


Originating  from  Accelerate  ter  Originating  fro* 

Scale  Factor  and  Mlsallgnaent  Errors  Accelerometer  Bias  Errors 


Originating  from  Accaleraaeter  Seal*  Factor  and 
Mlsallgnaent  Errors 


TABLE  A-2. 


INTEGRAL  FUNCTIONS  FOR  VELOCITY  ERROR  PROPAGATION  (concluded) 


Acceptable  Ratio 

of  Glabal  Rates 

JcosaS1n2sdt 

All  n 

J S1naCos2adt 

All  n 

-J  C0SaC0S2#dt 

All  n 

Jsin2oS1nedt 

n  *  2 

Jsin2oCcssdt 

All  n 

JsinaS1n2sdt 

All  n 

J'cosaSInZsdt 

All  n 

Jsin2«S1n2»dt 

No  n 

J S1n2aCos2sdt 

No  n 

Jcos2aS1n2adt 

No  n 

J Cos2aCos20dt 

No  n 

Jsin2aS1n2edt 

All  n 

Jcos2aS1n2Mt 

All  n 

Jsin2aS1n28dt 

All  n 

Jsin2aCos28dt 

All  n 

J  Sin  2«S1n2Mt 

n  t  1 

APPENDIX  B 


NAVIGATION  ERROR  MODEL 


The  navigation  error  eqnations  used  at  the  basis  for  performance  analysis  are 
defined  in  this  appendix,  together  with  a  discretization  procedure  for 
representing  the  error  model  in  the  form  of  difference  equations. 


The  navigation  error  model  assumed  is  a  '  ’vp-angle '  '  formulation,  which  is  a 
classic  set  of  navigation  error  equations  applied  to  either  strapdown  or 
platform  systems.  The  ' '^i-angle' *  error  model  may  be  expressed  in  the 
following  form: 


¥  =  ay  ( p  +  fl)  x  V 


(B-l) 


V  -  ea  -  ?  x  AL  -(20  *  b)  x  6Y  -  u2g  6R 


6V 


+  3<^(8R*S/R)E/R  +  5* 


&R  =  &V  -  *  x  5R 
where : 


(B-2) 

(B-3 


¥  *  attitude  error  vector 
&Y  3  velocity  error  vector 
6R  »  position  error  vector 

p.  »  angular  velocity  vector  of  local  vertical  reference  frame  relative  to 
an  earth-fixed  frame 

0  >  angular  velocity  vector  of  the  earth  relative  to  an  inertial  frame 
A^  a  vehicle  nongravitational  acceleration 
R  »  radius  vector  from  earth's  center  to  vehicle 
R  a  magnitude  of  £ 
ws  *  Schuler  frequency  ■  J  g/R 
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£.g  ”  gyro  output  error  vector 
el  =  accelerometer  output  error  vector 
*  gravity  error 

and  it  is  assumed  that  all  vectors  are  expressed  with  components  in  a 
local-vertical  navigational  reference  frame. 

The  total  gyro  and  accelerometer  error  vectors  are  expressed  in  terms  of 
individual  sensor  errors  as 


(HXg 

(B-4) 

OCX, 

(B-5) 

where: 

Xg  »  vector  of  gyro  errors  (bias  plus  random  noise) 

X,  *  vector  of  accelerometer  errors  (bias  plus  random  noise) 

C  =  transformation  matrix  relating  inertial  sensor  assembly 
reference  frame  to  the  local-vertical  frame 
G  =  coefficient  matrix  that  transforms  the  collection  of  accelerometer 

errors  into  a  net  acceleration  error  in  the  sensor  assembly  reference 
frame 

H  =  coefficient  matrix  that  transforms  the  collection  of  gyro  errors 

into  a  net  angular  rate  error  in  the  sensor  assembly  reference  frame 
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If  the  inertial  navigator  is  strapdown,  then  C  is  defined  to  be  the  direction 
cosine  matrix  relating  the  vehicle  body  axes  to  the  local-vertical  frame.  If 
the  inertial  navigator  is  of  the  platform  type,  the  matrix  C  is  the  identity 
matrix.  Thus,  Equations  B-l  through  B-5  are  applicable  to  both  types  of 
navigation  system,  with  the  appropriate  choice  for  the  matrix  C.  The  inertial 
navigation  system  error  equations  given  by  Equations  B-l  through  B-5  can  be 
conveniently  expressed  in  the  generalized  partitioned  form: 


Y 

X 


r. 

F„  . 

Y 

G.  „ 

0 

ii 

12 

+ 

11 

1 1 

L° 

X 

0 

°22 

*2 

(B-6) 


where : 

Y  =  vector  of  navigation  errors 
X  =  vector  of  sensor  errors 

It  is  desirable  to  discretize  Equation  B-6  in  the  partitioned  form: 

Y  =4>,,Y  +4>,_X  +r 

~n  lln-1  12  n-1  *  j 


^n  ^22-n-l  +  £  2 


where : 

(1  and  (2  are  appropriate  noise  vectors. 

The  discretization  of  the  sensor  states  normally  presents  no  problem  since  it 
is  commonly  assumed  that  sensor  errors  are  adequately  represented  by  either  of 
the  following  two  forms 

* 

b  =  t|  -  “  (first-order  Gauss-Markov  process) 
b  =  r)  (random  walk) 
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where: 


b  =  sensor  bias 

v  -  correlation  tine  of  process 
i)  =  white  noise 

If  the  sensor  error  is  a  first-order  Ganss-Markov  process,  the  discrete 
transition  and  noise  variance  parameters  (<b,D  are  given  by: 

<t>=  e-t/x 

02(1_e-2t/T) 

where  a  is  the  steady-state  standard  deviation,  and  x  is  the  transition  time. 
When  the  sensor  error  follows  a  random  walk  type  variation,  the  parameters  are 
<|>*  l 

r=  Nt 


where  N  is  the  spectral  intensity  of  the  white  noise  function,  i)  and  is 
expressed  in  convenient  units  as 


tandard  deviation  of  chan« 
60 


The  navigation  error  transition  matrix  can  be  discretized  by  using  a  Piccard- 
type  expansion  and  then  solving  by  numerical  integration.  The  result  is  given 
in  Tables  B-l  and  B-2  for  a  comprehensive  set  of  sensor  errors.  If,  as  in  this 
case,  the  gyro  and  accelerometer  misalignments  (non-orthogonalities)  are  given 
as  constant  random  errors,  only  nine  misalignment  angles  need  be  defined. 

Either  the  gyro  triad  or  the  accelerometer  triad  may  be  designated  as 


TABLE  B-l.  SELF-DEPENDENT  NAVIGATION  STATE  TRANSITION  MATRIX 
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the  ''reference''  triad  and  assigned  three  misalignment  angles,  with  six 
misalignment  angles  for  the  remaining  triad.  In  this  study,  the  gyros  have  six 
misalignment  angles,  and  the  accelerometers  have  three. 

The  discretization  algorithms  in  Tables  B-l  and  B-2  are  valid  over  times  that 
are  small  compared  to  the  Schuler  period  and  the  sensor  correlation  times.  The 
transition  and  noise  matrixes  can  be  extended  over  any  time  by  the  chain  rules: 

4>(t  ,  t  )  =  <t>  <t>(t  ,  t  ) 

n  o  n  n-1  o 

T 

Ht  ,  t  )  =  d>  r(t  .  ,t  )d>  +r 

n  o  Y j  o  td  n 

where : 

<p n  =  incremental  transition  matrix  for  interval  tn_i  to  tn 

fn  =  incremental  noise  covariance  matrix  for  interval  tn_j  to  tn 

The  covariance  matrix,  P,  of  the  error  states  can  be  reconstructed  at  any  time 
from : 

Pn  =  ^(tn*  to)po<t>T  <tn«to)  +  T  ( tn,  t0) 

The  error  model  defined  by  Equations  B-l  through  B-5  also  applies,  with  minor 
modification,  to  alignment  and  calibration.  The  principal  difference  is  that 
6R  is  a  fixed  error  and  has  no  variation  with  time.  Also,  since  the  error 
equations  presume  small  angles,  a  coarse  level  and  azimuth  alignment  is  implied 
as  a  prerequisite  for  their  use. 


In  addition  to  the  sensor  biases,  the  navigation  error  equations  are  driven  by 
two  additional  sources  of  randcn  error:  1)  sensor  random  noise,  and  2)  gravity 
errors.  These  additional  error  sources  are  accounted  for  by  means  of  the 
incremental  noise  covariance  matrix  defined  in  Table  B-3.  The  white  noise 
originating  from  the  gyros,  as  characterized  by  the  spectral  intensity  matrix 
Ng,  leads  to  the  attitude  random  walk.  The  white  noise  originating  from  the 
accelerometers  (absent  in  the  present  application)  is  characterized  by  the 
spectral  intensity  matrix  Nt  in  Table  B-3. 


i 


101 


iSt 


.ii: 


•tWi  Vi  aA  aLAiAaAsJmil  •h'l  %  'i  •TiJ 


TABLE  B-3 .  NAVIGATION  STATES  NOISE  OOVAKIANCE  MATRIX 


* 

SV 

6R 

<0 

qu  *  /CNgCTdt 

q12  *  /qllAT<lt 

"13  ‘  r"l2dt 

<5V 

q21  “  q12 

q22  *  ^Aq12dt 
+  (/Aql2dt)T 

+  (N#  +  Ny)t 

q23  ~ 

+  /q22dt 

5R 

q31  3  q13 

q32  “  q23 

"33  *  ;"23dt 
♦  (/q23dt)T 

Ng  =  Spectral  intensity  matrix  of  gyro  white  noise 


N  *  Spectral  intensity  matrix  of  accelerometer  white  noise 
3 


Ny  »  Spectral  intensity  matrix  of  white  noise  originating  in  local- 
vertical  frame  due  to  gravity  errors 


C  *  Sensor  to  local  vertical  transformation  matrix 


-a. 


a 

|-a 


‘1*  V 


h 


local-vertical  components  of  nongravitational  acceleration 
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APPENDIX  C 


SENSOR  QUANTIZATION  ERROR  MODEL 


The  discrete  output  of  inertial  instruments  causes  an  error  due  to 
quantization.  Tha t  is,  if  the  true  integral  of  the  sensed  rate  or 

acceleration  over  a  given  iteration  interval  is  not  equal  to  an  integer  number 
of  pulses,  the  fractional  pulse  will  be  "saved''  and  output  over  the  next 
interval.  No  fractional  pulses  are  ever  lost  and  the  output  of  the  instruments 
could  be  accumulated  indefinitely,  with  a  maximum  error  of  one  pulse.  This 
property  povides  the  basis  for  mathematically  defining  quantization  error. 

Define  the  quantization  error  over  the  nth  iteration  interval  as  qn.  Also, 
define  the  accumulated  quantization  error  as  qn.  Then: 

nn  =  nn-l  +  <ln  (C-l) 

The  cumulative  quantization  error,  qn,  constitutes  a  random  uncorrelated 
sequence  with  uniform  distribution.  Rearranging  Equation  C-l,  the  quantization 
error  over  the  nth  iteration  interval  is  expressed  as 


qn  =  nn  ~  nn-i  (c-2) 

Thus,  quantization  error  is  characterized  as  the  difference  between  two 
successive  members  of  a  random  independent  set. 

ATTITUDE  ERROR  DUE  TO  GYRO  QUANTIZATION 

Consider  attitude  error  in  a  rotational  environment  with  the  gyros  skewed.  The 
attitude  matrix  differential  equation  for  this  case  is 


C  -  C  (Am) 


(C-3) 


where: 


C  *  orthogonal  direction  cosine  aatrix 

A  *  transf oxaation  aatrix  frca  skewed  sensor  axes  to  a  set  of  orthogonal 
sensor  reference  axes 
a  =  angular  rate  aeasnxed  by  gyros 

and  {v}  denotes  a  skew- syaae tr i c  aatrix  foraed  froa  the  coaponents  of  a 
vector  v. 


The  discrete  algoritha  for  updating  Equation  C-3  is 


C 

n 


C  ,  I  +  CA©|  + 


^(A6}2+... 


where: 


(C-4) 


(C-5) 


Directly  froa  Equation  C-4,  the  error  in  C  is 
6C  -  6C 


n  'n-1 
+  C 


'n-1 


£l  +  (A©)  +  ^{A0}2+  ,..J 
j^A60  +  ^{AB0){A0}  +  *{A0}(A50]  +  ...J 


(C-«) 


where  is  the  error  in  the  increaental  gyro  output  defined  by  Equation  C-5. 


Since  gyro  quantization  error  can  only  cause  a  rotational  error,  the  erroneous 
C  aatrix  can  be  expressed  as 


A 

C 


[» -  <v]< 
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where: 
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Equation  C-ll  defines  a  generalized  recursive  relationship  for  attitude  error 
propagation  which,  when  69  represents  gyro  quantization  error,  takes  the 
specific  fora: 


-  V  'ST,-!1 


(C-12 


Equation  (C-12)  can  be  expanded  as 


V  «  C  Atk 
“n  1  -*1 


+C2A 


aj-UCj-Cjl+C^Aiij 

*c3a%-[(c3-E2)+c^]*% 


+. . . 
+c 


Av[("Cn-Vl)+Vl]AVl 


which,  after  cancellations,  reduces  to: 


n-1 


S.  ’  C»AS.  '  El  ACi-l*!!i-l 


(c-i3: 


where: 


ACi_i  =  Cj  -  C*_i 


The  covariance  of  ^  is  determined  directly  from  Equation  C-13  as 

E(¥  VT)  -  C  AATC  V  +  2  AC.AAT  AC.V 
m  n  q  g  i«l  i  i  g 


(C-141 


where: 


E(n  nT) 


Equation  C-14  indicates  that  gyro  quantization  produces  a  bounded  attitude 
error  (first  term),  plus  an  unbounded  error  component  (second  term).  The 
unbounded  coaponent  will  grow  at  a  rate  proportional  to  the  rotational  velocity 
of  the  gyro  triad. 

Another  fora  for  expressing  attitude  error  growth,  which  is  aore  convenient  for 
covariance  analysis,  is  given  by  the  equation  set: 


V  *  C  Afl  +  u _ , 

— n  n  n  — n~l 

u  *  u  ,  -AC  An 
“n  " n-1  n  Ti 


(C-15) 


which  is  siaply  a  rearrangeaent  of  Equation  C-13  with  the  sumaation  expressed 
recursively. 

VELOCITY  ERROR  DOE  TO  GYRO  AND  ACCELEROMETER  QUANTIZATION 


The  computed  direction  cosine  aatrix,  C,  is  used  to  transfora  accelerations 
measured  in  the  acceleroaeter  axes  into  components  in  the  reference  frame: 

V  =  Ca  (C-16) 

where : 

V  -  velocity  in  reference  fraae 

a  »  measured  acceleration  in  sensor  fraae. 

It  is  assumed  that  the  accelerometers  fora  an  orthogonal  triad  that  coincides 
with  the  orthogonal  gyro  reference  axes. 


A  discrete  algorithm  for  updating  Equation  C-16  is 


(C-17) 


where  vn  is  the  incremental  accelerometer  output,  and  C  is  the  direction  cosine 
matrix  at  the  midpoint  of  the  computational  interval.  The  error  incurred  in 
accumulating  velocity  is 


5V 


6V  4  +  5C  v  +  C  6v 
— n-1  n-n  n  ~t 


6V  -  (l)Cv  +  C6v 
n-1  "Ti  nm  n  i 


-  ‘Vi  -  x  4  +  c»Sl» 

-  ‘Vi *  t *  V*. 


( C— 1 8) 


where: 

vr  *  C  v 
— n  nn 


and  where  the  misalignment  vector  V  is  defined  as 

?  -  l/2()?  +  *  .) 

~n  a  ti-  x 


which,  for  all  practical  purposes,  is  accurately  approximated  by  '? n  ■ 


Equation  C-18  can  be  decomposed  into  two  parts  (to  separate  the  contributions 
due  to  gyro  and  accelerometer  quantization  errors)  as  follows. 


<  -  <-i  * 

<  -  *. 


(C-19) 

(C-20) 
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Thus,  the  effect  of  accelerometer  quantization  on  velocity  error  hat  exactly 
the  same  form  as  the  effect  of  gyro  quantization  error  on  attitude  error.  The 
only  difference  is  that,  because  the  gyros  are  taken  to  be  skewed,  the  attitude 
error  equation  involves  the  transformation  matrix  A.  By  omitting  A  and  making 
the  substitution  (6V*-*¥)  Equations  C-13,  C-14,  and  C-15  are  immediately 
applicable  to  accelerometer  quantization  error  propagation. 

The  complete  model  of  gyro  and  accelerometer  quantization  effects  is  summarized 
in  Figure  C-l .  The  model  is  recursive  and  operates  at  the  attitude  and 
velocity  update  frequency. 

COVARIANCE  ANALYSIS  MODEL  FOR  PROPAGATION  OF  QUANTIZATION  ERRORS 


For  covariance  analysis,  we  need  to  express  the  effect  of  sensor  quantization 
error  over  periods  of  time  much  longer  than  the  iteration  interval  of  the 
update  equations.  The  covariance  matrix  can  be  determined  by  solving  the 
discrete  covariance  equation 


ATTITUDE  AND  VELOCITY  UPOATE  INTERVAL 
UNCORRELATED  NOISE  SEOUENCES 


Figure  C-l.  Deterministic  Model  of  Quantization- Induced  Error 


for  the  error  aodel  defined  by  Figure  C-l .  The  evaluation  of  Equation  C-21  is 
straightforward  and  results  in  the  solution  given  in  Figure  C-2.  For 
convenience,  the  snauaations  can  be  replaced  by  continuous  integrations, 
assuming  that  the  variations  in  v£  are  relatively  smooth.  Then,  we  can  write: 

n-1  t-T  T 

E  AC,  AAtAC_  =  T  /  CAA  C  dt 
k=l  k  K 

t-T 

=  T  Iq  C{Au}AAT{Am}TCTdt 

where  T  is  the  update  interval. 

The  covariance  solution  given  by  Figure  C-2  defines  the  variances  and  cross 
covariances  of  the  errors  due  to  sensor  quantization  over  any  length  of  time. 

In  practice,  the  solution  is  carried  out  only  over  periods  of  time  equal  to  the 
discretization  interval  of  the  error  equations.  The  noise  covariance  matrices 
E(V  6V^) ,  and  E(6V8VT)  then  provide  the  statistical  description  of  a 
discrete  uncorrelated  input  to  the  navigation  error  equations  over  the 
interval.  After  each  successive  discretization  interval,  the  solution  inherent 
in  Figure  C-2  is  started  anew. 
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ovariance  Model  for  Quantization-Induced  Error 


